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ABSTRACT 

Infrared Dark Clonds (IRDCs) harbor the earliest phases of massive star for¬ 
mation, and many of the compact cores in IRDCs, traced by millimeter contin¬ 
uum or by molecular emission in high critical density lines, host massive young 
stellar objects (YSOs). We used the Robert C. Byrd Green Bank Telescope 
(GBT) and the Karl G. Jansky Very Large Array (VLA) to map NH 3 and CCS 
in nine IRDCs to reveal the temperature, density, and velocity structures and 
explore chemical evolution in the dense (> 10 22 cm -2 ) gas. Ammonia is an excel¬ 
lent molecular tracer for these cold, dense environments. The internal structure 
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and kinematics of the IRDCs include velocity gradients, filaments, and possibly 
colliding clumps that elucidate the formation process of these structures and their 
YSOs. We find a wide variety of substructure including filaments and globules 
at distinct velocities, sometimes overlapping at sites of ongoing star formation. 
It appears that these IRDCs are still being assembled from molecular gas clumps 
even as star formation has already begun, and at least three of them appear 
consistent with the morphology of “hub-filament structures” discussed in the lit¬ 
erature. Furthermore, we find that these clumps are typically near equipartition 
between gravitational and kinetic energies, so these structures may survive for 
multiple free-fall times. 

Keywords: molecular data - ISM: clouds - (ISM:) dust, extinction - ISM: 
molecules - Stars: formation - radio lines: ISM 


Introduction 


Infrared Dark Clouds (IRDCs) are dense (> 10 5 cm -3 ) and cold (< 20 K) collections 
of dust and molecular gas, typically arranged in filamentary and/or globular structures with 
compact cores. IRDCs were first observed as regions of high extinction in silhouette against 
the galactic background infrared (IR) emission by the Infrared Astronomical Satellite (IRAS) 
(Wood et al. 1994), the Midcourse Space Experiment ( MSX ) (Carey et al. |1998; Egan et al. 


1998 Simon et al. 2006a), and the Infrared Space Observatory (ISO) (Hennebelle et al. 


2001). The ability to identify and study these objects has been greatly improved by the 


Spitzer Space Telescope, and in particular surveys with two of its instruments: the Infrared 


Array Camera (IRAC) (Fazio et al. 2004) and Multiband Imaging Photometer for Spitzer 


(MIPS) (Ri eke et al. 2004). 


The Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) (Benjamin 


et al. 2003 Churchwell et al. 2009) covered the region |6| < 1° and 10° < \t\ < 65° in all the 


IRAC bands (3.6, 4.5, 5.8, and 8.0 /irn) with 1"5 to 1"9 resolution, while the MIPS Galactic 
Plane Survey (MIPSGAL) ( |Carey et al. 2009) was a complementary survey in the 24 and 
70 /irn MIPS wavebands with 6" and 18" resolution, respectively. These surveys observed 
much of the galactic plane, imaged the structure of IRDCs at higher resolution than before, 
and revealed embedded young stellar objects (YSOs), typically 1 to 10 per IRDC. (Note 
that in this paper we use the general ternPYSO,” encompassing both protostars and H II 
regions.) An extensive catalog of Spitzer IRDCs is given by Peretto & Fuller (2009), of which 
80% were previously not identified from MSX images. Furthermore, Simon et al. (2006b) 
matched several hundred IRDCs to molecular clouds seen in 13 CO .7=1-0 by the Boston 
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University Galactic Ring Survey (BU-GRS) (Jackson et al. 2006), and thus determined 
velocities, kinematic distances, and physical properties of the population of clouds. The 
darkest clouds have very high column densities, as high as approximately 10 2i -10 25 


cm 


-2 


The terms “core” and 


“clump” are frequently used in the literature, but the meanings 
are not standardized. We use the term “core” to refer to an unresolved or marginally resolved 
overdense structure <0.1 pc across and tens of solar masses. We then use the term “clump” 
to refer to a resolved structure (e.g. projected area greater than three times the area of 


the beam) within an IRDC and is assigned by a clump deconvolution algorithm (see 14.1 
for description of algorithms used). “Average” results of spectral line fitting presented in 
this paper are averaged over clumps. We also occasionally refer to the velocity components 
within IRDCs when the IRDC’s emission primarily occurs at two or more distinct velocities 
with little emission at the intermediate velocities; a velocity component may have one or 
more clump associated with it. 

IRDCs are an active subject of study because they are nurseries of massive star formation 
and contain complex chemistry. It is now widely speculated that these objects contain the 
earliest stages of the formation of massive stars. Recent studies have revealed massive (>10 


Mq) YSOs ( 

Rathborne et al. 2005| Pillai et al. 2006 

Beuther & Steinacker 2007 

Wang et al. 

2008) and massive (10-1000 M 0 ) cores in IRDCs ( 

hathborne et al. 

2006, 

2007 

Swift 

2009 


by the presence of masers (indicating accretion disks or outflows), radio continuum emission 
(from ionized gas), or fitting spectral energy distributions (SEDs) consistent with massive 
YSOs across mid-infrared (MIR) and submillimeter wavelengths. Kim et al. (2010) found 


that about 13% of IRDC cores have YSOs identified by their IR or maser emission. Not all 
dense cores in IRDCs are obviously star-forming, which raises the question: “Are these cores 
going to form stars and are just too young, or is there something fundamentally different 
about these cores preventing star formation?” To answer this question, we must know about 
the physical properties of both starless and actively star-forming clumps in IRDCs. 


Rathborne et al. (2006) used the Institut de Radioastronomic Millimetrique (IRAM) 30 


m telescope to make 11" resolution 1.2 mm dust continuum maps of 38 IRDCs to investigate 
the structure and clumps of IRDCs traced by dust emission. They found that these clumps 
had similar masses, sizes, and densities as hot (> 50 K) clumps with massive YSOs, but were 
much colder (15-30 K), implying that the clumps were indeed an early evolutionary phase 
of massive star formation if the clumps were to subsequently collapse. The authors further 
suggested IRDCs may be precursors to star clusters, as they have masses comparable to 
young clusters and contain several compact clumps. Finally, they also asserted that better 
resolution was needed to distinguish better individual cores and investigate fragmentation. 
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Disentangling the substructure of IRDCs is important to determining their formation 
and how it relates to subsequent star formation. Early studies of filamentary clouds with low 
spatial and velocity resolution appear to show that such structures have radial density profiles 


of the form p oc r 2 (see for example 

Alves et al. 1998 

Lada et al. 1999 

). This density profile 

is consistent with the presence of a toroidal magnetic field ( 

Fiege & Pudritz 2000 

), which can 

prevent the filaments from expanding ( 

Contreras et al. 2013 

). Higher resolution studies show 


that IRDCs typically have complex substructure, spatially and kinematically. Some studies 
observe distinct velocity components within IRDCs, seen in shock tracers such as SiO, that 
may be due to protostellar outflows or cloud-cloud collisions (Sanhueza et ah 2013). Other 


studies observe collections of filaments with coherent velocity gradients consistent with gas 


flowing toward the most massive cores, perhaps guided by magnetic fields (Peretto et al. 


2014). Myers (2009) and hi et ah (2013) describe a “hub-filament structure,” in which gas 


flows along filaments to a central hub where star formation is ongoing. Many studies have 
already described more generally how the collision and accumulation of smaller clouds and 
gas flows can assemble molecular cloud complexes (see for example Blitz & Shu [1980; Dobbs 


2008 

Furukawa et al. 

2009 

Torii et al. 2011 

et al. 

2014; 

Fukui et al. 

2015 

)• 


Ammonia is an excellent probe of the molecular gas in IRDCs. It is a relatively abundant 
species, typically about 10 -9 to 10~ 7 fractional abundance relative to molecular hydrogen 
(Ragan et ah 2011 Chira et ah 2013). Unlike carbon-bearing species, NH 3 suffers minimal 


freeze out in these cold environments (Bergin & Tafalla 2007). The ratio of the (1,1) to (2,2) 
line strengths leads to a rotation temperature that has been shown to be a good tracer of 


the kinetic temperature of the gas over the typical range for IRDCs (Ho & Townes 1983). 
Furthermore, the critical density of NH 3 is about 10 4 cm -3 , well matched to the gas being 


studied (Tafalla et ah 2002). (Simon et ah (2006b) found typical average number densities 

to 10 4 cm -3 by observing 13 CO J=l-0 with 


for the majority of IRDCs to be about 10 3 cm 


a beam a few times larger than most of the IR extinction features, so when accounting for 
beam dilution we expect the dense gas associated with these features to be near or above the 
critical density for ammonia.) Ammonia has hyperfine structure, allowing determination of 
the optical depth in the NH 3 lines by comparing the observed hyperfine component strengths 
to those set by molecular physics (Ho & Townes 1983). The direct measurement of the 


temperature and optical depth allows for an unambiguous determination of the column 
density, perfect for environments in which the column can vary significantly from optically 
thin to optically thick in the same IRDC. 


Devine et al. (2011) observed the NH 3 (1,1) and (2,2) and the CCS (2 i-1q) transitions 


toward G19.30+0.07, and measured temperatures around 15-20 K, NH 3 column densities 
around 10 15 -10 16 cm -2 , and linewidths of about 2 km s -1 . They also observed that NH 3 
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and CCS generally had spatial distributions that were anticorrelated. Furthermore, they 
found that G19.30+0.07 was composed of four distinct clumps in three distinct velocity 
components. These clumps had masses of tens to hundreds of solar masses and were virially 
unstable against self-gravitating collapse. This complex substructure has also been investi¬ 


gated in other star-forming regions with physical properties similar to IRDCs. Hacar et ah 


(2013) additionally found that molecular filaments in Taurus could be deconvolved into a hi¬ 


erarchy of several sub-filaments at distinct velocities leading to core formation. (The smallest 
scale of the hierarchy, the “filament bundles,” have spatial widths and linewidths comparable 
to the separation of the bundles, about 0.1 pc and 0.5 km s^ 1 , so our study does not have 
the resolution to identify similar bundles in these IRDCs if they occur on the same scales as 
in Taurus.) Rosolowsky et al. (2008) also saw multiple velocity components in NH 3 toward 
dense cores in Perseus at 0.04 pc, 0.024 km s^ 1 resolution. These individual components 
often had such narrow linewidths (~0.1-0.3 km s -1 ) that turbulence was unlikely to con¬ 


tribute significantly to slowing collapse. Lu et al. (2014) studied 62 high-mass star-forming 


regions with Very Large Array (VLA) observations of the NH 3 (1,1) and (2,2) transitions. 
They found that parsec scale filaments were ubiquitous and often contained regularly spaced 
dense cores, though we note this regular spacing could be a result of using interferromet- 
ric observations. Furthermore, they suggested that the filaments could be supported by 
turbulence and found that the dense cores were near virial equilibrium. 


Ragan et al. (2011, 2012) performed a study of NH 3 (1,1) and (2,2) with the Robert 


C. Byrd Green Bank Telescope (GBT) and the Very Large Array (VLA) in six IRDCs at 
3'!7-8'!3 and 0.6 km s -1 resolution. They found that the majority of the gas had kinetic 
temperatures 8-13 K, indicating that protostcllar heating was not significant for most of 
the IRDC. Furthermore, they found that velocity fields were generally coherent across the 
IRDCs, with local (~0.1 pc) disruptions of a few km s -1 coincident with sites of local 
star formation. They argued that neither turbulence nor thermal pressure were sufficient 
to support the IRDCs and that the observed velocity structures were a result of ongoing 
collapse, fragmentation, and protostellar feedback. 


Ragan et al. (2013) further investigated the differences in filamentary and globular 


IRDCs and their hierarchical structure. They studied 11 IRDCs, covering a range of mor¬ 
phology and star formation activity, with Herschel and the Submillimetre APEX Bolometer 
Camera (SABOCA) instrument on the Atacama Pathfinder Experiment (APEX) 12 m tele¬ 
scope at 350 /irn, resolving structure down to ~0.1 pc. They performed a dendrogram 
analysis on the APEX data found that filamentary IRDCs tended to be more massive and 
have more hierarchical structure than clumpy IRDCs. This suggests that IRDCs may be 


divided into two relatively distinct morphological families. Ragan et al. (2014) identified 7 


Giant Molecular Filaments (GMFs) in which the molecular gas extend for ~100 pc including 
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IRDCs, IR bright structures, and more diffuse gas presumably enveloping these denser re¬ 
gions. The existence and structure of GMFs suggests that hierarchical structure may extend 
to even these large size scales. 


Single-dish data are necessary for determining whether sharp transitions at the edges of 
the dense structures are real and not a consequence of interferometers filtering the extended 
emission. Fortunately, the NH 3 (1,1) and (2,2) and the CCS (2!-l 0 ) transitions may be 
observed simultaneously with the GBT. It is also important to remember that we likely do not 
resolve the smallest substructure in these IRDCs and that interferometric data are even more 
important. The combination of interferometric and single-dish radio data from the VLA and 
the GBT gives us the high resolution we need to separate better individual cores and other 
structures compared with previous work, while also not resolving out emission. The limited 
number of studies in this regime suggest that the substructure is present because IRDCs 
represent a phase in which star formation proceeds even as the IRDC itself is still being 


constructed from colliding clumps and/or gas flowing into the IRDC along filaments (Myers 


2009 Li et al. 2013 Sanhueza et al.|2013[). A homogeneous study of several representative 


IRDCs combining total power and high resolution is thus necessary. For the current study, 
we mapped the NH 3 ( 1 , 1 ) and ( 2 , 2 ) and the CCS (2i-lo) transitions across nine IRDCs 
containing both ongoing star formation and more quiescent environments. 


We describe the IRDC sample, with distance determination and previous studies of these 
IRDCs, in (|2j Observations and data reduction are described in §[3j Our methodology is 
presented in (|4j including NH 3 spectral line fitting, clump deconvolution, and IR extinction. 
Results of investigating the kinematics, the spectral line fitting, and clump properties are in 
£j5j A discussion of the results relating to kinematics, structure, gravitational stability, and 
chemical evolution is presented in §[ 6 j and conclusions are in i|T| 


2. Sample Selection and Distance Determination 


This sample of nine IRDCs was originally selected in 2005 from extinction features in 
the 8 pm GLIMPSE images and chosen to cover a wide range of physical parameters. The 


IRDCs subsequently appeared in a catalog compiled by Simon et al. (2006a), identified by 


their high contrast against background emission in the 8.3 pm MSX images. The sample 
size is large enough to contain tens of cores identified by dust thermal continuum, including 
about half that are coincident with IR point sources indicating ongoing star formation, but 
is also small enough for each object to be studied in detail. The larger IRDCs show both 
quiescent, starless cores and YSOs in close proximity to each other. 
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We include both filamentary and globular morphologies, apparently starless and pro- 
tostcllar cores, and ranges of physical sizes, masses, linewidths, distances, and IR contrast. 


A comparison of this sample with all the IRDCs analyzed by Simon et ah (2006b) in the 


BU-GRS data is shown in Figure [I] (except G010.74-00.13, which is not covered by the BU- 
GRS). The BU-GRS is a large scale survey of the 110.2 GHz 13 CO J=l-0 transition in the 
disk of the Milky Way using the Five College Radio Astronomy Observatory (FCRAO) 14 
m single dish telescope. The publicly available data cubes have velocity resolution of 0.2 km 
s -1 , angular resolution of 46", and typical antenna temperature RMS sensitivity of 0.13 K 
(Jackson et al. 2006). This sample tends toward higher mass and column density the general 
population of IRDCs in the BU-GRS. 


We calculated the kinematic distances using the galactic rotation curve of Reid et ah 


(2009). They adopted a galactocentric radius R 0 = 8.4+0.6 kpc and a circular rotation speed 
©o = 254 ± 16 km s -1 kpc -1 , based on the results of their measured trigonometric parallaxes 
of massive star-forming regions. All of the regions in our sample were in the \t\ < 90° regime, 
so there was naturally a near-far distance ambiguity. Since IRDCs are seen in extinction, we 
assumed for the rest of this work that all of the regions were at the near kinematic distance, 


except for G034.43+00.24. Kurayama et al. (2011) performed parallax measurements of an 


HoO maser associated with G034.43+00.24 using Very Long Baseline Interferometry (VLBI) 
as part of the VLBI Exploration of Radio Astrometry (VERA) project and determined the 
distance to be 1.56lo ii kpc, less than half of the near kinematic distance of 3.55lg 


+0.36 tj. 
36- iL 


should be noted that the distance to this IRDC is still under debate. Foster et al. (2012) used 


two near-infrared (NIR) extinction methods on two different NIR data sets to derive four 
extinction distances, all of them consistent with the near kinematic distance and three of 
them inconsistent with the maser distance (the distance estimate with the largest uncertainty 
was consistent with the maser distance). They further raise questions about the size of the 
uncertainties on the maser distance given the quality of the model fit to the data and the 
low declination of the source. We adopt the maser distance in this work, noting that sizes 
and masses reported for this IRDC may be factors of several times larger if the kinematic 
and extinction distances more accurate. A summary of the sample is presented in Table [TJ 


3. Observations and Data 
3.1. GBT Observations 

Data from the GBT provide total flux on large spatial scales, which is not measured 
by our VLA observations. Dates of GBT observations for Project IDs AGBT05C_014 and 
AGBT12B_283 are listed for each source in Table [2] The data are dual polarization taken 
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Fig. 1.— A comparison of physical parameters of the IRDCs in this sample to all of those 


in the BU-GRS (Jackson et al. 2006) as analyzed by Simon et al. (2006b). The black points 


represent IRDCs from this sample (except G010.74-00.13, which was not covered by the 
BU-GRS). 


Table 1. IRDC Sample 


IRDC 

R.A. (J2000) 

Deck (J2000) 

UjSR 

D 


hh:mm:ss.s 

ddmmrss 

(km s' 1 ) 

(kpc) 

G010.74-00.13 

18:09:45.9 

-19:42:04 

29.0 

3 4fi+ 0 - 46 

G022.56-00.20 

18:32:59.6 

-09:20:08 

76.8 

4 55 +a26 

G024.60+00.08 

18:35:39.7 

-07:18:52 

53.4 

o A c+0.32 
,3 - 4O -0.34 

G028.23-00.19 

18:43:30.2 

-04:12:56 

80.0 

4 50 +0 ' 30 

^• du -0.30 

G031.97+00.07 

18:49:30.6 

-00:48:18 

96.0 

5 60 +a42 

O.OU— 0 42 

G032.70-00.30 

18:52:06.2 

-00:20:26 

90.2 

C 00+0.45 
' J - ZO -0.39 

G034.43+00.24 

18:53:18.7 

01:25:51 

58.0 

1 56+ 012 

G035.39-00.33 

18:57:09.4 

02:07:48 

44.5 

9 ns+0- 37 

z - yo -0.38 

G038.95-00.47 

19:04:08.4 

05:09:12 

42.2 

9 on+0.40 
^•ot | -0.40 
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in frequency switching mode using a 5 MHz shift. Observations in 2005 are in PointMap 
mode with the single pixel K-band receiver, while data in 2012 are in on-the-fly (OTF) 
mapping mode using beams 1 and 4 of the 7 possible beams in the K-band Focal Plane 
Array (KFPA). All observing configurations include simultaneous observation of the NH 3 
( 1 , 1 ) and ( 2 , 2 ) inversions lines and the CCS (2i-lo) line, with rest frequencies of 23.6945 
GHz, 23.72263 GHz, and 22.34403 GHz, respectively. The expected beam FWHM at these 
frequencies is approximately 32 / . / Hourly pointing scans indicated the actual FWHM of the 
beam during observing varied between 32" and 34", and the pointing correction was typically 
6 " 


The spectral resolution of the raw data varied between observations with different spec¬ 
tral setups, so all data were smoothed to the limiting velocity resolution of approximately 
0.15 km s' 1 . System temperatures typically varied between 40 and 60 K. Absolute flux den¬ 
sity calibration of the 2012 observations was tied to beam nodding scans of Venus on 2012 
December 05. We utilized the Green Bank Telescope Interactive Data Language (GBTIDL) 
procedures venusmodel, venuscal, and venuscalget, documented in GBT memo ^275^ 
to model the antenna temperature of Venus at a given date, time, and elevation, and to 
calculate the aperture efficiency. This model corrected for the distance to Venus and any 
beam dilution if it was unresolved, but did not account for variations across the surface of 
Venus. The determined aperture efficiency, tja, of 0.64 was in good agreement with the value 
of 0.67 from the Ruze equation and the GBT sensitivity calculator for the NH 3 (1,1) line, 
although we note both of these are approximations for extended sources and there will be 
some variation at different frequencies. From this value, we also calculated the main beam 
efficiency to be r ] m b = 1.37 tja — 0.88 as noted in “Calibration of GBT Spectral Line Data in 
GBTIDL v2.1.”0 

The primary uncertainty in the efficiency estimates come from the possibility of pointing 
errors in the observations of Venus, which will tend to cause us to underestimate the aperture 
efficiency. Given typical pointing scan corrections, we do not expect this effect to be bigger 
than about 30%, which would result in the flux densities in our final maps being 30% higher. 
The good agreement between the measured aperture efficiency and the prediction form the 
Ruze equation implies that the errors are not this large. 

Inspection of the observations taken in 2012 shows that while the relative flux density 
scales of dates 2012 November 14 and 2012 December 5 agree, observations on day 2012 
November 5 have a lower flux density scale. This effect is seen in the maps of the same regions 


'https://safe.nrao.edu/wiki/bin/view/GB/Knowledge/GBTMemos 
"http://www.gb.nrao.edu/GBT/DA/gbtidl/gbtidl_calibration.pdf 



Table 2. GBT Observations 


Date 

Session # 

IRDC 

Zenith Opacity 

Project ID AGBT05C_014 

2005 Oct 19 

1 

G010.74-00.13 

0.123 

2005 Oct 19 

2 

G038.95-00.47 

0.115 

2005 Oct 23 

3 

G010.74-00.13 

0.128 

2005 Oct 30 

6 

G022.56-00.20 

0.057 

2005 Oct 31 

7 

G022.56-00.206 

0.083 

2005 Oct 31 

7 

G024.60+00.08 

0.084 

2005 Oct 31 

7 

G032.70-00.30 

0.084 

2005 Nov 02 

8 

G032.70-00.30 

0.067 

2005 Nov 03 

8 

G038.95-00.47 

0.067 

2005 Nov 17 

9 

G032.70-00.30 

0.039 

2005 Nov 17 

9 

G038.95-00.47 

0.04 

2005 Nov 18 

11 

G028.23-00.19 

0.058 

2005 Nov 19 

11 

G038.95-00.47 

0.056 

2005 Dec 24 

12 

G028.23-00.19 

0.066 

Project ID AGBT12B_283 

2012 Nov 05 

2 

G034.43+00.24 

0.048 

2012 Nov 14 

3 

G034.43+00.24 

0.043 

2012 Dec 05 

4 

G031.97+00.07 

0.041 

2012 Dec 05 

4 

G035.39-00.33 

0.040 



on different days, the peak antenna temperature in the Focus scans of our pointing source 
1851+0035 (observed approximately every hour), and calibration observations of DR21. To 
determine the relative flux density scale, we averaged the results of the fractional difference in 
the maps and the Focus scans. We excluded DR21 in this calculation because the measured 
ratio may be affected by pointing errors since we neither mapped DR21 nor performed a 
pointing correction on it. We measured the flux density scale on 2012 November 5 to be 0.7 
relative to that of the other days, so we adopted values of the aperture efficiency and main 
beam efficiency of 0.46, and 0.63, respectively, to bring data from that day into agreement 
with the others. This discrepancy may be a result of atmospheric effects, dish surface effects, 
or pointing or focus errors not accurately corrected by the standard calibration procedure. 

No flux calibrators were observed during the 2005 observations, so we adopted 0.58 and 
0.8 as the aperture efficiency and main beam efficiency, respectively, because the aperture 
efficiency at these frequencies in 2005 was typically 91% of what is is now (T. Hunter 2015, 
private communication). We inspected the peak antenna temperature of the focus scans on 
the pointing sources across multiple epochs to look for variations in the relative flux density 
scales, assuming the flux densities of these sources to be stable over the approximately two 
month timescale of these observations. The most frequently observed pointing source was 
1733-1304, being observed every day of the observations in 2005 except 2005-12-24, when only 
1743-0350 was used. Additionally, 1751+0939 was observed sporadically on 2005 October 
19, 2005 November 17, and 2005 November 18. Assuming these calibrators were constant in 
flux density, we found that the aperture efficiency only changed significantly over the course 
of 2005 December 24, possibly affecting the flux density scale by about 20% across the maps 
of G028.23-00.19, though this may also be due to the intrinsic variability of the calibrators. 

Initial calibration of the frequency switching data and conversion to main beam tem¬ 
perature was performed with the GBT pipeline for the KFPA provided by the NRAO. The 
pipeline queried the weather data from the GBT weather forecast database to determine the 
zenith opacity for each observing block. The pipeline also accepted values for the aperture 
efficiency and the main beam efficiency as inputs, determined as above. 

After the initial calibration, we used GBTIDL to fit and subtract baselines from the 
spectra. Many of the baselines contained complex shapes with large amplitudes. We averaged 
the spectra for each individual combination of scan number, IF, polarization, and beam, then 
fit the line-free regions with a high order polynomial. The baseline shapes varied considerably 
between different IFs, beams, and polarizations, and showed moderate variation with time 
among different scans. The baselines were typically approximately linear in the vicinity of the 
line, though not always. We determined that the variation between individual integrations 
was within the noise of an individual integration, so that we could average over whole scans. 
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For the 2012 observations taken in OTF mode, a single scan consisted of one row or column 
in the map. Since the 2005 observations were taken in PointMap mode, each scan was a 
single pointing. Baselines in the 2005 observations were generally less complex, so we used 
lower order polynomials at this stage. 

Imaging was performed also using the GBT pipeline for the KFPA, which uses Parsel- 
Tongue to call AIPS tasks from Python. Data were gridded using the SDGRD task. Since 
lower amplitude baseline effects were noticed on shorter timescales, we additionally per¬ 
formed lower order polynomial fitting and subtraction on an individual line of sight basis 
in the final imaged data cubes using the imcontsub task in the Common Astronomy Soft¬ 
ware Applications (CASA) package (http://casa.nrao.edu). By combining two baseline 
subtraction methods, we could address both the time variability in neighboring scans before 
the imaging step while also taking care to get the flattest baselines possible in the final 
data cubes. The data were converted from main beam temperature to Jy beam -1 by mul¬ 
tiplying by (2/cB^mb)/(A p ?/ J 4) = 0.483, where ks is the Boltzmann constant and A p is the 
physical collecting area of the GBT. Typical RMS noise in the final data cubes was 25-40 
mJy beam -1 per 0.15 km s -1 channel. We estimate our absolute flux density uncertainty to 
be approximately 30%. 


3.2. VLA Observations 

NH 3 (1,1) and (2,2) and CCS (2i-l 0 ) emission in our sample of IRDCs was observed with 
the VLA in D configuration, primarily in 2005 and 2006 (NRAO Proposal ID AD0516), with 
observations of CCS in G028.83-00.19, G031.97+00.07, and G034.43+00.24 in 2007 (NRAO 
Proposal ID AD0556). The 3.125 MHz bandwidth covered only the main and innermost 
satellite hyperhne lines of the NH 3 (1,1) and (2,2) transitions. The spectral setup had 
24.414 kHz resolution. Different numbers of pointings per IRDC were used to fully cover the 
highest opacity regions with the approximately 1'9 (FWHM) primary beam, as determined 
from the 8 /jrn images. G031.97+00.07 and G034.43+00.24 each were observed with five 
pointings (though only four were observed in CCS), G028.83-00.19 was observed with two 
pointings, and the remaining IRDCs were observed with one pointing each. J1820-2528 was 
used as the phase and amplitude calibrator for G010.74-00.13, while J1832-1035 was used 
for G022.56-00.20 and G024.60+00.08, and J1851+0035 was used for the remaining IRDCs. 
3C286 (J1331+3030) and 3C48 (J0137+3309) were used as flux density calibrators. Table [3] 
summarizes the parameters of the VLA observations, including bandpass calibrators, beam 
sizes, and flux-density-to-temperature conversion factors. 

All data were calibrated using AIPS. The 2007 CCS observations were conducted during 


Table 3. VLA Observations 


IRDC 

Project 

ID 

Bandpass 

Calibrator 3, 

Synthesized Beam 
Major Axis x Minor Axis 
(") x (") (pc)x(pc) 

PA 

(deg) 

Flux to Tb 

(K Jy 1 ) 

NH 3 (1,1) 23.6945 GHz 

G010.74-00.13 

AD0516 

J1851+0035 

4.8x3.4 

0.081x0.057 

-4.6 

134.6 

G022.56-00.20 

AD0516 

J1832-1035 

4.3x3.3 

0.095x0.073 

-2.3 

155.3 

G024.60+00.08 

AD0516 

J1832-1035 

4.2x3.3 

0.070x0.055 

7.3 

158.2 

G028.23-00.19 

AD0516 

J1851+0035 

4.0x3.6 

0.087x0.079 

-9.6 

154.2 

G031.97+00.07 

AD0516 

J1851+0035 

3.9x3.6 

0.106x0.010 

30.3 

153.6 

G032.70-00.30 

AD0516 

3C273, J1851+0035 

3.8x3.8 

0.096x0.096 

64.9 

149.3 

G034.43+00.24 

AD0516 

3C273, J1851+0035 

3.8x3.7 

0.029x0.028 

79.2 

156.3 

G035.39-00.33 

AD0516 

3C273, J1851+0035 

4.2x3.6 

0.061x0.052 

62.6 

143.8 

G038.95-00.47 

AD0516 

3C273, J1851+0035 

3.9x3.6 

0.053x0.049 

87.0 

154.2 

NH 3 (2,2) 23.72263 GHz 

G010.74-00.13 

AD0516 

J1832-1035 

4.7x3.4 

0.079x0.057 

-7.6 

134.0 

G022.56-00.20 

AD0516 

3C273, J1832-1035 

4.4x3.3 

0.097x0.073 

-6.3 

149.3 

G024.60+00.08 

AD0516 

3C273, J1832-1035 

4.2x3.3 

0.070x0.055 

1.7 

154.8 

G028.23-00.19 

AD0516 

3C273, J1851+0035 

4.2x3.5 

0.092x0.076 

-23.1 

149.3 

G031.97+00.07 

AD0516 

3C273, J1851+0035 

3.9x3.6 

0.106x0.098 

-9.1 

154.3 

G032.70-00.30 

AD0516 

3C273, J1851+0035 

4.1x3.7 

0.104x0.094 

-45.2 

144.3 

G034.43+00.24 

AD0516 

3C273, J1851+0035 

3.8x3.7 

0.029x0.028 

-69.3 

153.0 

G035.39-00.33 

AD0516 

3C273, J1851+0035 

4.2x3.7 

0.061x0.053 

64.9 

141.0 

G038.95-00.47 

AD0516 

3C273, J1851+0035 

4.0x3.7 

0.054x0.050 

-78.4 

148.9 

CCS (2i-l 0 ) 22.34403 GHz 

G010.74-00.13 

AD0516 

3C454.3 

5.4x3.2 

0.091x0.054 

15.8 

140.3 

G022.56-00.20 

AD0516 

3C454.3 

4.4x3.5 

0.097x0.077 

8.0 

160.9 

G024.60+00.08 

AD0516 

3C273, 3C345, 3C454.3 

4.3x3.6 

0.072x0.060 

-0.5 

158.5 

G028.23-00.19 

AD0556 

3C345 

4.8x3.5 

0.105x0.076 

-5.3 

147.5 

G031.97+00.07 

AD0556 

3C345 

5.0x3.4 

0.136x0.092 

-9.2 

142.5 

G032.70-00.30 

AD0516 

3C273, 3C454.3 

4.3x3.6 

0.109x0.091 

-46.2 

159.8 

G034.43+00.24 

AD0556 

3C345 

4.8x3.4 

0.036x0.026 

-10.8 

147.1 

G035.39-00.33 

AD0516 

3C84, 3C345, 3C454.3 

4.1x3.6 

0.059x0.052 

-49.2 

162.5 

G038.95-00.47 

AD0516 

3C84, 3C345, 3C454.3 

4.1x3.6 

0.056x0.049 

-55.1 

166.7 


a Bandpass calibrators were selected by the best bandpass solution for each observing date. CCS obser¬ 
vations often also used the phase and amplitude and/or flux density calibrators to improve the bandpass 
solution. 3C273 is also known as J1229+020, 3C454.3 is also known as J2253+161, 3C345 is also known 
as J1642+398, and 3C84 is also known as J0319+041. 
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the EVLA upgrade, and eight of the twenty-seven antennas had been converted to EVLA 
antennas. Data were obtained using both the VLA and EVLA antennas. Doppler tracking 
could not be used during the upgrade period, so the AIPS task CVEL was applied during 
calibration to correct for motion of the Earth relative to the LSR during the observations. 
In the 22 GHz data sets, the uncertainty of the absolute flux densities was estimated to be 
10%. We estimated the absolute positional accuracy to be better than 1" 

Before imaging, data were Hanning smoothed, bringing the effective velocity resolution 
to approximately 0.6 km s -1 . Continuum subtraction was performed with the CASA task 
uvcontsub on all the data for G034.43+00.24 and the CCS observations of G031.97+00.07, 
as there were significant continuum point sources seen in these four dirty images. The data 
were imaged in CASA using the clean task in mosaic mode with natural weighting of the 
visibilities, deconvolved with 0.33 km s -1 channels and 077 pixels. We also employed the 
multiscale capability of clean to include clean components approximately one and three 
times the size of the synthesized beam, as well as the standard point-like components. The 
CASA task pbcor was used to apply a primary beam correction. The mosaics were imaged 
to the 35% power level relative to the peak of the mosaic primary beam response. 

The synthesized beam in the cubes is 3"-5". The RMS noise is 1-4.5 mJy beam -1 per 
0.33 km s -1 channel. The conversion between from flux density in Jy beam -1 to brightness 
temperature, T B (K), for the various data cubes are listed in Table [3j 


3.3. Combining Single-dish and Interferometric Data 

Successful combination of single dish and interferometric data relies upon good astromet¬ 
ric alignment and calibration between the two data sets. The estimated positional accuracy 
of our GBT observations, 6", while small compared to the GBT beam, were larger than the 
VLA beam. This allowed for positional errors that might introduce significant artifacts in 
our combined maps. For example, using the GBT as a clean model for the VLA with a 
positional offset between the two can cause negative features neighboring the emission in 
the combined image. To address this issue, we computed the amount of positional shift 
required to align the GBT and VLA data by smoothing our VLA-only (1,1) cubes to the 
GBT angular resolution and smoothed the GBT-only (1,1) cubes to the VLA velocity resolu¬ 
tion. We then determined the positional shift required to align the maxima in the emission. 
Shifts were restricted to integer numbers of pixels (6" each) in the native GBT cubes. The 
largest total shift was 18", or approximately half a GBT beam, for G035.39-00.33. All other 
IRDCs required shifts less than 10", and three IRDCs (G032.70-00.30, G034.43+00.24, and 
G038.95-00.47) did not require any shift. 
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The GBT data were combined with the VLA data in a two-stage process. First, the 
GBT cube was used as a starting model for the VLA data in the CASA task clean, and then 
cleaning proceeded as for the VLA only, described above in 1 3.2 The resulting cube was 
then combined with the GBT cube again, using the CASA task feather. Before feathering, 
the VLA primary beam response was applied to GBT cube with the same cutoff as clean to 
ensure that emission in the GBT cube that was outside the VLA footprint was not included in 
the Fourier transform. Primary beam correction was then applied to the resulting feathered 
cubes. No relative calibration factor was applied during feathering, implicitly assuming that 
there was no systematic offset in the absolute calibration of the datasets. 


Feathering Fourier transforms the two images, adds the two filtered Fourier cubes, and 
trans- forms back. Feathering two images is a simple way to combine the two data sets, but 
can be sensitive to the shape of the tapering function and to relative calibration uncertainties. 
Using a cleaned image that has the single dish data as a starting model, rather than the 
cleaned VLA image alone, mitigates these effects because the cleaning process can correct any 
potential issues with the feathered image that might be inconsistent with the interferometric 
visibilities. Combining the data sets by using a clean model alone, however, is also sensitive 
to calibration errors and can result in a cube with more total flux than the single dish alone, 
which is not physical. Feathering these resulting cubes with the single-dish brings the total 
flux in the final combined images back into agreement with the single-dish data. 


Using both methods produced a result that is as consistent as possible with the both 
the GBT and VLA data taken individually. The RMS noise was approximately the same as 
the VLA-only cubes. The synthesized beams, and thus the conversion between brightness 
temperature and flux density, were exactly the same as for the VLA data alone (Table [3]). 

As a check on the flux in our combined cubes, we smoothed them to the resolution of the 
GBT and divided by the GBT cubes. This should have produced a cube of values close to 1, 
as the GBT is sensitive to total flux. We found that within regions of significant emission, 
our smoothed, combined cubes were almost always within the 20% flux uncertainty of the 
GBT. Our final cubes were therefore within the uncertainty of recovering the correct total 
flux. Comparing to the VLA data cubes alone, the morphology of the combined data cubes 
is similar on size scales approximately the size of the synthesized beam, however using the 
VLA alone misses most the extended, diffuse emission. The fluxes of clump-sized sources 
seen in the VLA are tens of percent lower than what is seen in the combined images because 
of the contribution of this diffuse gas across the IRDCs. The total flux in the VLA cubes 
is typically less than half of the total flux measured by the GBT, and thus the combined 
images. An image showing the comparison between GBT-only, VLA-only, and GBT+VLA 
images is shown in Figure [2} 



16 


(a) GBT (b) VLA (c) GBT+VLA 



40s 36s 32s 28s 24s 18h49m20s 

Right Ascension (J2000) 


Fig. 2.— A comparison of the NH 3 (1,1) data for the GBT-only, VLA-only, and combined 
images. Images of G031.97+00.07 are integrated from 94 km s ^ 1 to 97 km s _1 . The VLA 
footprint is plotted on the GBT image as a reference. The VLA and combined images have 
both been corrected for the response of the primary beam and are shown on the same flux 
scale (all three flux scales are linear). Aside from the improved resolution over the GBT 
along, the combined image also recovers more extended emission than the VLA alone. The 
combined image recovers the total flux of the GBT image, and recovers more flux than the 
VLA for even relatively compact, marginally resolved sources. 
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Methods 


4.1. Clump Deconvolution 

Detailed knowledge of the kinematic and spatial structure of these IRDCs provides 
constraints on the NH 3 spectral line fitting (see 14.3). It is readily apparent in the data 


that regions containing at least two strong, distinct velocity components, if not three or 
four, are common. Attempts to fit a single velocity component to the NH 3 spectra in 
these regions result in poor fits with unphysically large linewidths and optical depths, so 
multiple components must be included. By performing clump deconvolution on the data 
before fitting, it is trivial to determine the number of components to include at each line of 
sight. Sophisticated clump deconvolution algorithms make use of the data cube and its noise 
properties as a whole, and so are more robust determinations of the number of components 
to fit than any method making this determination on a line of sight (pixel-by-pixel) basis. 
Additionally, identifying significant, coherent emission in the data allows us to make very 
good initial estimates of the central velocities and velocity widths at each line of sight for 
each component by making first- and second-order moment maps (velocity field and velocity 
dispersion, respectively) restricted to emission within the identified clump. Finally, using 
these clumps provides a straightforward and physically motivated way to analyze the physical 
parameter results from the fitting for different substructures within the IRDCs. 

We perform clump deconvolution on the NH 3 (1,1) data using the cprops package 


described by Rosolowsky & Leroy (2006). We use the main hyperfine component of the 


( 1 , 1 ) line because it typically has the highest signal-to-noise ratio and the ( 2 , 2 ) line does 
not trace the coldest, and often lowest column density, gas. The velocity offsets between 
cospatial velocity components are sometimes comparable to the NH 3 hyperfine splitting, so 
the deconvolution must be performed carefully. We use the GBT and VLA combined data 
cubes before primary beam correction so the noise is roughly constant across the images 
to generate the clump assignment cubes. The VLA primary beam response was applied to 
the GBT data before combination with the uncorrected VLA data specifically to allow the 
cubes to be used for clump deconvolution in this way; applying the VLA primary beam 
response to the final combined data cubes would be equivalent. The initial mask only 
includes voxels (single elements in the position-position-velocity cubes, akin to pixels in 
position-position images) with values greater than 7<r, and then the mask is expanded to 
include voxels above 5cr that are connected to the initial mask via a path only passing through 
significant emission. The cprops algorithm then identifies “kernels” in the data, consisting 
of local maxima significantly above a “merge level” - the contour value at which multiple 
kernels are connected in the data and the significant emission cannot be uniquely assigned 
to one kernel over another. The cprops documentation refers to this unassigned emission 
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as the “watershed.” We further restrict the list of kernels to those that lead to clumps with 
projected area on the sky greater than three times the synthesized beam area, i.e. only 
well-resolved clumps are assigned as independent structures so that we can accurately probe 
their physical properties. Kernels must also result in clumps that extend over at least 3 
velocity channels (recall channels are 0.33 km s -1 ). Kernels that lead to these unresolved 
clumps are rejected before the final cprops assignments are determined, so the voxels, and 
thus emission, from these clumps is free to be reassigned to another clump or the watershed. 
This rejection ensures that all clumps are resolved, effectively rejecting cores. This algorithm 
produces clump assignments from our data that map excellently to clumps discerned by eye. 


An alternate method in the cprops package uses the clumpfind algorithm (Williams 


et al. 1994) to assign emission to the same kernels list as the standard cprops algorithm, 


but proceeds by assigning voxels to clumps in discrete contours continuing all the way to the 
noise floor. Assignment degeneracies in this method are broken by evaluating the proximity of 
voxels to the clump peaks in position-position-velocity space. This method has the advantage 
of assigning all of the significant emission into clumps, however the assignment cubes tend 
to have a “patchwork” appearance that likely does not represent physically accurate clump 
boundaries. The clumpfind algorithm is also fairly sensitive to the input parameters that 
determine the contouring scheme used to make the clump assignments. This will in turn 
affect the clump parameters we calculate in £j6| The cprops algorithm is generally less 
sensitive than clumpfind to the inputs. The specific assignments can be altered by varying 
the parameters, but the typical clump sizes, aspect ratios, etc. are not significantly affected. 


For our study, we start with the standard cprops assignments, and then further used 
the clumpfind assignments to assign the watershed emission. The result is an assignment 
cube in which all of the significant emission is assigned to exactly one clump, and the 
patchwork assignments from the clumpfind algorithm are restricted to the weakest emission. 
The clump-averaged properties discussed in this work will be dominated by the strongest 
emission, and thus by the cprops assignments. Results of the deconvolution are presented 


in £5.1 and discussed in £6.4 


4.2. IR Extinction 

IRDCs are initially identified by their apparent extinction in MIR images, so it is natural 
to use their contrast with the background emission as a measure of physical properties in 
the IRDC. In the simplest terms, a greater contrast, i.e. darker cloud compared to the 
background, indicates a greater column of dust. It is rather straightforward to generate 
extinction maps if one has an estimate of the background, and with a few assumptions the 
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dust mass surface density may also be calculated across a map. 

Any method of estimating this optical depth must account for the following complica¬ 
tions: (1) we do not have a direct measurement of the background IR emission at the location 
of the IRDC, so it must be estimated from an irregular, varying background measured off 
the IRDC while attempting to avoid contamination from foreground objects, such as other 
dark clouds and bright nebulae; (2) there exists foreground emission from the dust between 
the IRDC and the observer; and (3) IRDCs often contain embedded point sources that do 
not probe the full column of the IRDCs and contaminate our contrast measurements. 

We adopt a modified version of the Large Median Filter (LMF) method presented by 
to map the 24 /iin optical depth, r 2 4 (Um , using the Spitzer MIPSGAL 
data. This method is summarized in detail in Appendix [A] In our optical depth maps, the 
IRDCs typically peak at optical depth of about 0.25-0.5. G028.23-00.19 has noticeably higher 
optical depth than the other IRDCs by about a factor of two, and peaks at an approximate 
optical depth of 1. Further results of this analysis are presented in §5.2[ 


Butler & Tan (2009 


4.3. Ammonia Spectral Line Fitting 


We fit the spectra along individual lines of sight in our combined GBT and VLA data. 
The (1,1) and (2,2) lines were fit simultaneously. The fitting routine was written in Python 
using the nrapfit package, which performed a least-squares (Levenberg-Marquardt algo¬ 
rithm) fit and returned the fit parameters and the full covariance matrix. We include the 
full hyperhne structure of NH 3 from the components and intrinsic strengths documented by 


Kukolich (1967). The limited bandwidth of the VLA data restricted our fit in practice to the 


main and inner satellites components of the spectral lines. We simultaneously fit the central 
velocity (u C! lsr); velocity FWHM (An), total optical depth in the (1,1) component (r 0 (l, 1), 
abbreviated henceforth as r 0 ), excitation temperature (T ex ), and rotation temperature (T R ). 
Note that here we use Tq for the total opacity in the entire (1,1) line, which is precisely a 
factor of two larger than the total opacity in the main component of the (1,1) line. 

There is a natural degeneracy between the excitation temperature and optical depth 
if they are low. However, the main hyperhne component is typically optically thick in this 


study (see 15.3), and so the degeneracy is broken. Since we observe lines that are so optically 


thick (to > 5), our fitted values of To may only be lower limits. The fitting routine is capable 
of handling multiple velocity components simultaneously and independently along the same 
line of sight (see Section 4.1). 


The details of the spectral line fitting routine and the determination of physical pa- 
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rameters, including the kinetic temperature, T K , and the column density of ammonia are 
presented in Appendix [B] A list of key parameters for the analysis, including the parameters 
used in the spectral line fitting, is given in Table |4} Results of the spectral line fitting are 
presented in §5.3[ An example spectrum and fit of a single line of sight with two distinct 
velocity components is shown in Figure [3j 



a(J2000) = 18 h 49 m 33 s , S{ J2000) = -00 ° 47'48 


v lsr ( km s 1 ) 


Fig. 3.— An example spectrum of the NH 3 (1,1) and (2,2) lines with fit to the spectrum. 
This spectrum is from one line of sight toward G031.97+00.07 with two distinct velocity 
components. The data are the dashed lines and the fit is the solid line. The (2,2) spectrum 
and fit have been vertically offset by -3 K for clarity. The brightness temperature scale 
assumes a beam filling factor of unity. 
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Table 4. Selected Physical Parameters of Interest 

Parameter Description 

Line of Sight Parameters 

Calculated Directly from the Radio Data: 

Moment 0 (integrated intensity) of the NH 3 (1,1) line 
(i ; ( 1 , 1 )) Moment 1 (intensity-weighted velocity field) of the NH 3 (1,1) line 
(^,( 1 , 1 )) Moment 2 (intensity-weighted velocity dispersion) of the NH 3 (1,1) line 
Ices Moment 0 (integrated intensity) of the CCS line 

Calculated Directly from the IR Data: 

T 9 A/i .m Optical depth at 24 /irn from IR extinction 

From NH 3 (1,1) and (2,2) Spectral Line Fitting 8, : 

u c Central Doppler shifted NH 3 frequency 
Av NH 3 frequency FWHM 
r 0 Total optical depth in the NH 3 (1,1) line 
T ex NH 3 excitation temperature 
Tr NH 3 rotation temperature 
X 2 Goodness of fit statistic 

Calculated from Spectral Line Fitting Results 8 : 

Tk Kinetic Temperature 
Vc,lsr Central NH 3 LSR velocity 
Av NH 3 velocity FWHM 
dobs Observed line of sight velocity dispersion 

< 7 ii n e Line of sight velocity dispersion corrected for spectral resolution 
a t NH 3 Thermal velocity dispersion 
unt NH 3 Nonthermal velocity dispersion 
A4 Nonthermal (turbulent) Mach number 
N( NH 3 ) Total column density of NH 3 
1V(H 2 ) Column density of molecular Hydrogen 
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5. Results 

Images of the sample from the Spitzer Space Telescope and the Herschel Space Obser¬ 
vatory with the GBT and VLA data are shown in Figures [4} [5j [6j [7J [8} [9j [10| [TTJ and 
The NH 3 (1,1) distribution generally traces the IR extinction very closely and peaks around 
infrared point sources and, to a lesser extent, the IR extinction peaks. The NH 3 (2,2) dis¬ 
tribution is more compact and generally correlates with stronger (1,1) emission. The CCS 
emission had systematically lower signal-to-noise ratio (< 5<r) than the NH 3 and is typically 
only marginally detected. The CCS that is observed, however, does not typically cover the 
full spatial extent of the NH 3 . Extended CCS emission may be resolved out by the VLA, 
however the GBT images in Figures |4|fT2l show that frequently the strongest CCS emission 
is not cospatial with the strongest NH 3 emission. 
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5.1. Clump Deconvolution Results 

A summary of the clump properties is presented in Table [5] Coordinates and velocities 
are given at the peaks of the emission. The effective radius, _R e ff, is the radius of a circle on 
the sky that has the same projected area as the full extent of the clump using the clump 
assignment from cprops, which ranges from approximately 0.02 pc to 0.28 pc in our sample, 
with a median of about 0.16 pc. The physical resolution (i.e. synthesized beam) of the (1,1) 
cubes ranges from 0.02 pc to 0.1 pc, so clumps are typically well resolved. 
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Table 4—Continued 

Parameter Description 

Clump-by-clump Parameters 

Directly from cprops for each clump: 

R.A. (J2000) R. A. of the peak NH 3 (1,1) main beam temperature 
Deck (J2000) Deck of the peak NH 3 (1,1) main beam temperature 

-R c ff Effective radius of a circle with the same projected area as the clump 
A 0 (Initial) Aspect ratio projected onto the sky 

Calculated clump-by-clump: 

M c i Clump mass from total NH 3 and AT(NH 3 ) 
kf,s P h Spherical free-fall time 
fff,cyi Cylindrical free-fall time along the axis 
M v ir, sp h Spherical virial mass 
Af v ir, C yi Cylindrical virial mass 
«vir,s P h Spherical virial parameter 
«vir,cyi Cylindrical virial parameter 

B cr Minimum magnetic held strength to support clump against collapse 
Whole IRDC Parameters 
D Distance to IRDC 

/f or e Fraction of IR galactic emission that is foreground to the IRDC 
Mir Mass lower limit from 24 /jrn extinction 
M 13 co Mass lower limit from BU-GRS 13 C0 (,7=1-0) 

M„m Mass from BGPS 1.12 mm emission 

a These parameters may also be evaluated as clump-averaged values, weighted by the 
(xfed) 1//2 statistic over the clump. Both line of sight and clump-averaged values separate 
different velocity components that overlap spatially. 
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G010.74-00.13 



M I I I I I I 1 I I I I I I I I I I I I I I I I I I I I I 1 l I I I I I I I I I I I I I 


3.6\im 4.5pm 8.0pm 


8.0pm 24pm 70pm 


70pm 160pm 350pf 


(d) GBT NH,(1,1) 


(e) GBT NH,(2,2) 


(f) GBT CCS(2j—1 0 ) 


1 I 1 I I I I 1 1 I I 1 1 I 


(g) GBT+VLA NH,(1,1) 


(h) GBT+VLA NH,(2,2) 


(i) GBT+VLA CCS(2 1 -1 0 ) 


4.5pm 8.0pm 24pm 


Right Ascension (J2000) 


Fig. 4.— Spitzer Space Telescope and Herschel Space Observatory images of G010.74-00.13, 
along with our GBT and VLA data. The IRDC is seen against the diffuse mid-infrared 
emission of the galaxy toward the shorter wavelengths (panels (a) and (b)), though embed¬ 
ded infrared point sources are visible. At longer wavelengths (panel (c)), the IRDC becomes 
visible in the thermal emission of the dust. The highest 24 /jrn optical depth, r 2 4 /im , that we 
calculate in this IRDC is 0.42. A yellow scalebar representing 1' with the physical size at 
our adopted distance is shown at the lower right of (c). The GBT footprint is overplotted in 
dashed magenta in panels (a)-(c). Panels (d)-(f) show the NH 3 (1,1), NH 3 (2,2), and CCS 
emission, respectively, as seen by the GBT. An arbitrary linear color stretch is shown with 
white contours showing the signal-to-noise ratio at 5a, 10cr, 20cr, and 40c. The size of the 
GBT beam is shown in the lower left of each panel, and the VLA footprint is overplotted 
in dashed gray. Panels (g)-(i) show Spitzer images of G010.74-00.13 with contours of the 
combined GBT and VLA NH 3 (1,1), NH 3 (2,2), and CCS emission, respectively. The con¬ 
tours show the signal-to-noise ratio at 5a, 25a, and 50a. The NH 3 generally traces the MIR 
extinction and peaks near IR point sources. The NH 3 (2,2) is generally less extended than 
the (1,1), as it traces warmer gas. The CCS, though weak, typically follows the highest 
extinction parts of the IRDCs when it is detected. The combined footprint of the GBT and 
VLA observations (their intersection) is shown in dashed white. The synthesized beam is 
shown at the lower left of each panel. 










Declination (J2000) 
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G022.56-00.20 


03s 33m00s 57s 18h32m54s 4.5|im 8.0pm 24|im 

Right Ascension (J2000) 


Fig. 5.— Same as Figure [4] but for G022.56-00.20. The highest value of r 2 4 Atm is 0.27. 
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45s 42s 39s 18h35m36s 

Right Ascension (J2000) 


G024.60+00.08 


4.5|im 8.0pm 24pm 


Fig. 6.— Same as Figure [4] but for G024.60+00.08. The highest value of r 24jUm is 0.47. The 
green circle in panel (a) shows the location of the extended green object (EGO) identified 
by iCyganowski et ah (2008). 
















Right Ascension (J2000) 


Fig. 7.— Same as Figure [4] but for G028.23-00.19 with GBT+VLA contours at 5cr and 25cr. 
The bright point source is an unrelated, foreground late-type star (Bowers & Knapp 1989;). 
The highest value of r 2 4 Mm is 0.89. 
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G031.97+00.07 


42s 36s 30s 24s 18h49ml8s 4.5|im 8.0|im 24|im 

Right Ascension (J2000) 


Fig. 8.— Same as Figure [4] but for G031.97+00.07 with GBT contours at 5a, lOcr, 20cr, 40cr, 
80cr, and 160cr, and with GBT+VLA contours at 5a, 25a, 50a, and 100a. The highest value 


of T 2 4 Mm is 0.43. The black crosses in panel (c) are H 2 0 masers reported by Wang et al. 


(2006). The cyan circle shows the location of the H II region identified by Urquhart et al. 
(2009). 
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G032.70-00.30 


12s 10s 08s 06s 04s 02£8h52m00s 4.5|im 8.0pm 24pm 

Right Ascension (J2000) 


Fig. 9.— Same as Figure [4] but for G032.70-00.30 with GBT contours at 5cr and 10cr, and 
with GBT+VLA contours at 5a. The highest value of t?*,,™ is 0.36. 
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G034.43+00.24 


3.6pm 4.5pm 8.0pm 8.0pm 24pm 70pm 70pm 160pm 350pm 



28s 24s 20s 16s 18h53ml2s 4.5pm 8.0pm 24pm 

Right Ascension (J2000) 


Fig. 10.— Same as Figure [4] but for G034.43+00.24 with GBT contours at 5a, 10a, 20a, 
40a, 80a, and 160a, and GBT+VLA contours at 5a, 25a, 50a, and 100a. The highest value 
°f T 94 /Cm is 0.14. The green circles in panel (a) show the location of the extended green 


objects (EGOs) identified by Cyganowski et al. (2008). The black crosses in panel (c) are 

H 2 O masers reported by Wang et al. 

(2006). The cyan circle shows the location of the H II 

region identified by Urquhart et al. ( 

2009). 
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G035.39-00.33 


15s 12s 09s 06s 18h57m03s 4.5pm 8.0pm 24pm 

Right Ascension (J2000) 


Fig. 11.— Same as Figure [4] but for G035.39-00.33 with GBT contours at 5 a, 10cr, and 20cr. 
The highest value of T 24 Mm is 0.44. 
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G038.95-00.47 


12s 09s 06s 19h04m03s 4.5pm 8.0pm 24pm 

Right Ascension (J2000) 


Fig. 12.— Same as Figure[4]but for G038.95-00.47 with GBT+VLA contours at 5 a and 25cr. 
The highest value of T 24 /jm is 0.34. 











Table 5. Clump Properties 


Parent Peak Coordinates cprops Herschel Spectral Line Fitting Spitzer 


IRDC a 

R.A.(J2000) 
hh:mm:ss 

Decl. (J2000) 
dd:mm:ss 

^LSR 
(km s -1 ) 

Rett 

(pc) 

A 0 

70 fin i 
source? 

(^line) 

(km s -1 ) 

(to) 

(Tk) 

(K) 

(7V(NH 3 )) 

(10 23 cm“ 2 ) 

(t2Aiith) 

G10 

18:09:45 

-19:42:29 

28.50 

0.17 

2.1 

N 

1.0 

7.3 

14.1 

0.8 

0.35 

G10 

18:09:45 

-19:42:07 

28.50 

0.16 

1.2 

Y 

1.0 

8.3 

15.1 

1.1 

0.34 

G10 

18:09:44 

-19:42:06 

29.16 

0.11 

2.1 

N 

0.7 

7.9 

14.0 

0.6 

0.37 

G10 

18:09:46 

-19:41:47 

29.82 

0.16 

2.5 

Y 

1.0 

8.1 

14.7 

0.9 

0.27 

G22 

18:33:00 

-09:19:59 

74.84 

0.15 

3.0 

Y 

1.3 

6.7 

16.5 

1.0 

0.20 

G22 

18:33:01 

-09:19:56 

76.82 

0.09 

1.6 

N 

1.6 

6.0 

15.7 

0.9 

0.19 

G22 

18:33:00 

-09:20:04 

78.14 

0.09 

2.8 

N 

1.6 

5.1 

16.3 

0.7 

0.22 

G24 

18:35:40 

-07:18:36 

52.49 

0.12 

1.4 

Y 

1.4 

4.5 

19.8 

0.8 

0.19 

G24 

18:35:39 

-07:18:55 

53.15 

0.17 

1.5 

N 

0.8 

6.1 

14.0 

0.6 

0.28 

G24 

18:35:40 

-07:19:06 

53.81 

0.09 

1.2 

N 

0.6 

6.1 

14.2 

0.4 

0.23 

G24 

18:35:40 

-07:18:31 

53.81 

0.10 

1.2 

N 

1.0 

4.5 

18.1 

0.6 

0.15 

G24 

18:35:41 

-07:18:03 

53.81 

0.03 

3.4 

N 

0.3 

8.4 

11.7 

0.6 

0.06 

G24 

18:35:41 

-07:19:08 

54.14 

0.05 

1.6 

N 

0.4 

5.6 

16.0 

0.4 

0.19 

G28 

18:43:31 

-04:13:23 

79.15 

0.09 

1.5 

N 

1.2 

12.2 

13.9 

1.5 

0.84 

G28 

18:43:30 

-04:13:09 

79.81 

0.08 

1.3 

N 

1.0 

9.5 

13.6 

0.8 

0.77 

G28 

18:43:30 

-04:12:59 

80.14 

0.08 

1.8 

N 

1.5 

6.2 

12.3 

0.8 

0.70 

G28 

18:43:29 

-04:12:34 

80.14 

0.05 

1.8 

N 

0.9 

4.6 

13.8 

0.4 


G28 

18:43:29 

-04:12:27 

80.14 

0.06 

2.9 

N 

1.1 

3.5 

14.6 

0.4 

0.44 

G28 

18:43:30 

-04:12:30 

80.80 

0.06 

4.4 

N 

1.1 

9.9 

13.1 

0.9 

0.63 

G28 

18:43:30 

-04:13:32 

81.13 

0.05 

1.7 

N 

0.8 

8.5 

14.0 

0.7 

0.67 

G28 

18:43:29 

-04:12:51 

81.13 

0.04 

1.8 

N 

1.2 

13.0 

11.4 

1.2 

0.56 

G31 

18:49:29 

-00:48:53 

93.85 

0.18 

1.3 

N 

1.0 

7.6 

15.3 

0.9 

0.15 

G31 

18:49:26 

-00:49:52 

94.18 

0.09 

2.0 

N 

0.6 

7.3 

12.9 

0.5 

0.29 

G31 

18:49:31 

-00:47:09 

94.18 

0.10 

1.2 

N 

2.3 

4.3 

18.9 

0.9 

0.12 

G31 

18:49:28 

-00:48:30 

94.51 

0.17 

1.5 

N 

1.2 

4.7 

18.1 

0.6 

0.18 

G31 

18:49:29 

-00:48:03 

94.51 

0.17 

1.9 

Y 

1.3 

5.5 

18.1 

0.7 

0.18 

G31 

18:49:30 

-00:47:23 

94.51 

0.07 

1.7 

N 

0.7 

6.6 

16.4 

0.6 

0.17 

G31 

18:49:30 

-00:48:18 

94.84 

0.07 

2.9 

N 

1.0 

6.1 

16.5 

0.5 

0.28 

G31 

18:49:32 

-00:48:04 

94.84 

0.17 

1.5 

N 

1.7 

5.4 

17.1 

0.9 

0.29 

G31 

18:49:30 

-00:46:46 

94.84 

0.08 

1.4 

N 

1.6 

4.4 

16.9 

0.7 

0.14 

G31 

18:49:32 

-00:46:29 

94.84 

0.08 

2.3 

N 

1.8 

2.8 

19.5 

0.6 


G31 

18:49:26 

-00:50:28 

95.50 

0.06 

3.4 

N 

1.1 

6.8 

12.7 

0.8 

0.28 

G31 

18:49:31 

-00:47:17 

95.50 

0.10 

1.3 

N 

1.8 

3.6 

17.3 

0.6 

0.25 

G31 

18:49:33 

-00:47:15 

95.50 

0.16 

1.7 

N 

2.2 

5.4 

18.4 

1.1 

0.14 

G31 

18:49:24 

-00:50:07 

95.83 

0.09 

3.3 

N 

1.0 

5.6 

13.8 

0.5 

0.29 

G31 

18:49:34 

-00:46:44 

95.83 

0.28 

1.1 

Y 

1.8 

5.2 

18.9 

1.2 

0.01 

G31 

18:49:28 

-00:49:39 

96.49 

0.18 

2.0 

N 

0.9 

7.2 

13.2 

0.8 

0.28 

G31 

18:49:33 

-00:47:33 

96.49 

0.25 

2.1 

N 

1.7 

6.6 

16.7 

1.2 

0.25 

G31 

18:49:25 

-00:50:09 

96.82 

0.08 

4.7 

N 

1.5 

7.4 

12.2 

1.0 

0.25 

G31 

18:49:31 

-00:46:46 

98.47 

0.06 

2.1 

N 

1.3 

4.1 

18.0 

0.5 


G31 

18:49:32 

-00:46:59 

99.46 

0.18 

1.3 

Y 

2.0 

6.3 

18.0 

1.2 

0.12 

G31 

18:49:31 

00:46:41 

99.46 

0.12 

2.5 

N 

0.9 

5.4 

17.6 

0.6 

0.05 

G32 

18:52:03 

00:21:09 

89.83 

0.04 

1.3 

N 

0.5 

4.9 

10.5 

0.5 

0.31 

G32 

18:52:05 

00:20:54 

90.16 

0.08 

3.7 

Y 

0.8 

5.2 

14.3 

0.5 

0.25 
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Additionally, cprops performs an elliptical fit to each clump with a variable position 
angle to determine sizes along the major and minor axes. The aspect ratio, Aq, is calculated 
from the cprops results as the ratio of the major axis length to the minor axis length, which 
is confirmed visually to be an accurate method. The aspect ratios of most of the clumps are 
approximately 1-2, but extend to about 3 to 4 for the more filamentary clumps, while the 
highest aspect ratio in our sample is 6. 

The effective radii and aspect ratios we measure will be sensitive to the orientations of 
elongated structures with respect to the line of sight. The values we report for R c s and A 0 
will be lower limits, and so clumps will be systematically more filamentary than discussed 
here. We further discuss the projection effects in Appendix [C| 


5.2. IR Extinction Results 


The IR optical depth, T^m, is calculated as described in l ]4.2| and then averaged over the 
full extent of the clump on the sky. Since the extinction is calculated from two-dimensional 
images, we cannot accurately separate the extinction from multiple components along the 
same line of sight, however the majority of individual pixels in the extinction maps are 
included in at most one velocity component. Values of r 2 4 /tm range from approximately 0 to 
0.85. Clumps at the lower end of this distribution are those dominated by IR point sources 
or other emission, such that we cannot reliably determine an IR optical depth. Since none 
of the optical depths averaged over a whole clump are greater than 1, we can conclude that 
while the centers of clumps may be very optically thick, there are at least significant outer 
portions of the clumps that are susceptible to heating by external or internal IR radiation 
fields. They are, however, still opaque to optical and ultraviolet radiation since the extinction 


at 550 nm is approximately 50 times greater than at 24 /mi (Draine 2003) 


5.3. Spectral Line Fitting Results 

Plots of the results of pixel-by-pixel NH 3 spectral line fitting are shown in Figure [I3j as 
well as averages over clumps (also in Table [5]) and whole IRDCs. The parameter averages 
are weighted by the square root of the reduced chi-squared statistic from the fit, 

^ = Gw) ’ (i) 

where y 2 is the goodness of fit statistic and u-dof is the number of degrees of freedom, i.e. 
the difference between the number of data points and the number of fit parameters (five 








Tabic 5—Continued 


Parent 

IRDC a 

Peak Coordinates 
R.A. (J2000) Deck (J2000) 
hh:mm:ss dd:mm:ss 

^LSR 
(km s -1 ) 

cprops 

Reft A 0 
(pc) 

Herschel 
70 fim 
source? 

(^line) 

(km s _1 ) 

Spectral Line Fitting 

(to) (T k ) (JV(NH 3 )> 

1 (K) (10 23 cm" 2 ) 

Spitzer 

(T24 M m) b 

G32 

18:52:07 

00:20:06 

90.82 

0.05 

2.0 

N 

0.4 

3.6 

12.6 

0.3 

0.32 

G34 

18:53:17 

01:23:45 

55.86 

0.10 

2.1 

N 

1.4 

5.7 

19.7 

0.8 

0.01 

G34 

18:53:17 

01:24:32 

56.19 

0.08 

2.0 

N 

1.6 

5.3 

18.9 

1.3 


G34 

18:53:17 

01:26:52 

56.19 

0.05 

1.8 

N 

0.6 

6.5 

12.8 

0.5 

0.01 

G34 

18:53:18 

01:24:51 

56.52 

0.08 

1.1 

Y 

2.6 

4.4 

24.6 

2.0 


G34 

18:53:18 

01:26:37 

56.52 

0.04 

2.8 

N 

0.4 

6.6 

11.9 

0.5 

0.01 

G34 

18:53:16 

01:26:30 

56.85 

0.04 

2.6 

Y 

1.3 

6.5 

13.6 

0.9 


G34 

18:53:15 

01:26:58 

56.85 

0.02 

2.2 

N 

0.4 

9.4 

13.1 

0.6 


G34 

18:53:17 

01:23:59 

57.18 

0.02 

1.9 

N 

1.2 

4.7 

13.8 

0.7 


G34 

18:53:21 

01:23:08 

57.51 

0.02 

2.5 

N 

0.4 

8.6 

11.2 

0.6 

0.00 

G34 

18:53:18 

01:25:25 

57.51 

0.08 

1.4 

Y 

2.3 

4.3 

27.9 

1.8 


G34 

18:53:19 

01:26:27 

57.51 

0.09 

2.5 

N 

1.0 

5.6 

14.1 

0.7 

0.04 

G34 

18:53:19 

01:24:29 

57.84 

0.14 

2.7 

Y 

1.7 

6.1 

20.6 

1.6 

0.01 

G34 

18:53:18 

01:25:37 

58.17 

0.11 

1.2 

N 

1.5 

5.5 

19.8 

1.3 


G34 

18:53:19 

01:23:13 

58.50 

0.05 

2.3 

N 

1.0 

6.7 

14.7 

0.8 


G34 

18:53:18 

01:25:13 

58.50 

0.06 

1.6 

N 

1.5 

6.0 

21.3 

1.7 


G34 

18:53:19 

01:27:49 

58.50 

0.13 

2.9 

Y 

1.1 

6.3 

13.9 

0.9 

0.06 

G34 

18:53:20 

01:28:21 

59.49 

0.09 

1.4 

Y 

1.2 

5.2 

20.3 

1.1 

0.05 

G34 

18:53:21 

01:26:57 

59.82 

0.02 

1.8 

N 

0.9 

5.4 

13.8 

0.7 

0.01 

G35 

18:57:08 

02:08:25 

42.53 

0.09 

2.1 

N 

0.8 

7.1 

12.8 

0.7 

0.20 

G35 

18:57:10 

02:07:36 

44.84 

0.14 

1.6 

N 

0.8 

7.9 

14.2 

0.8 

0.34 

G35 

18:57:09 

02:07:08 

45.17 

0.04 

1.9 

N 

0.4 

6.7 

11.8 

0.5 

0.27 

G35 

18:57:09 

02:08:18 

45.17 

0.09 

1.4 

Y 

1.0 

5.2 

14.4 

0.6 

0.23 

G35 

18:57:07 

02:08:41 

45.17 

0.03 

2.0 

N 

0.7 

5.2 

14.0 

0.7 

0.25 

G35 

18:57:08 

02:07:52 

45.50 

0.10 

1.3 

Y 

0.9 

7.4 

13.5 

0.8 

0.32 

G35 

18:57:09 

02:07:52 

45.50 

0.11 

1.8 

Y 

1.0 

5.3 

14.5 

0.6 

0.26 

G35 

18:57:08 

02:08:03 

45.50 

0.13 

1.4 

Y 

0.7 

8.1 

13.9 

0.7 

0.29 

G38 

19:04:07 

05:08:47 

42.16 

0.10 

1.3 

N 

1.2 

5.2 

17.6 

0.9 

0.23 

G38 

19:04:10 

05:09:14 

42.16 

0.04 

2.0 

N 

0.6 

4.9 

22.9 

0.3 

0.22 

G38 

19:04:06 

05:08:50 

42.49 

0.07 

1.9 

Y 

1.3 

4.0 

13.2 

0.6 

0.16 

G38 

19:04:08 

05:09:30 

42.49 

0.10 

1.7 

Y 

1.3 

3.8 

15.1 

0.5 

0.16 

G38 

19:04:10 

05:08:43 

42.82 

0.04 

1.1 

N 

0.6 

4.5 

13.6 

0.5 

0.24 

G38 

19:04:10 

05:08:54 

42.82 

0.05 

1.6 

N 

1.1 

3.9 

15.0 

0.5 

0.10 

G38 

19:04:08 

05:08:56 

42.82 

0.10 

2.0 

Y 

1.3 

3.2 

15.3 

0.5 

0.26 

G38 

19:04:08 

05:09:09 

42.82 

0.11 

1.1 

N 

1.0 

3.4 

17.1 

0.4 

0.27 

G38 

19:04:07 

05:09:46 

42.82 

0.08 

1.4 

Y 

1.1 

4.9 

18.3 

0.8 

0.27 

G38 

19:04:10 

05:09:05 

43.15 

0.04 

1.3 

N 

0.8 

4.2 

15.0 

0.4 

0.17 


a Parent IRDC (values of / fore given in parentheses): G10 = G010.74-00.13 (0.095), G22 = G022.56-00.20 (0.173), G24 = G024.60+00.08 
(0.127), G28 = G028.23-00.19 (0.202), G31 = G031.97+00.07 (0.328), G32 = G032.70-00.30 (0.289), G34 = G034.43+00.24 (0.058), G35 
= G035.39-00.33 (0.131), G38 = G038.95-00.47 (0.131). 

b We cannot compute T 24 ^ m for some clumps because IR emission covers their entire angular extent or the contrast with the background 
is too low. 
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per velocity component). This is a natural choice for the weighting because the relative 
uncertainties of the fit parameters scale with (Xred) 1 ^ 2 - Experimentation shows that the 
averages are not very sensitive to the exact weighting scheme; for example, the results 
change by only a few percent if weighting by yf ed is used instead of (x^ed) 1 ^- 


The typical uncertainties for single lines of sight and single velocity components for 
u Ci lsr, An, r 0 , T ex , and T R are 0.04 km s -1 , 0.09 km s -1 , 1.4, 0.58 K, and 1.1 K, respectively. 
The velocity FWHM, kinetic temperature, optical depth in the (1,1) line, and column density 
are generally correlated. G034.43+00.24 (cyan) is notable as having the highest An, Tk, 
and 1 V(NH 3 ) in the sample, dominated by the emission surrounding the bright YSOs with 
outflows as traced by masers and extended green objects (EGOs). EGOs are bright, resolved 
4.5 //m sources often attributed to shocked molecular hydrogen emission in protostellar 
outflows (Cyganowski et al. 2008). G028.23-00.19 is also notable for being apparently starless 
throughout, and is colder and has a higher tq than the rest of the sample. 


The clump-averaged values quoted for the optical depth, r 0 , and the kinetic temperature 
of the gas, T Kl are averages weighted by the reduced (x^ed ) 1 ^ 2 va l ue from the spectral line 
fitting, and can be separated kinematically by individual clump even when they overlap 
spatially. Most of the clumps have an average NH 3 optical depth tq greater than 5 and all 
are greater than 1, indicating that these clumps are typically optically thick in NH 3 emission. 
We note that G028.83-00.19 shows significantly higher NH 3 optical depth than the rest of 
the sample, matching its high 1R optical depth. 


5.3.1. Linewidths 


In further analysis we adjust the values of the velocity dispersion to account for ob¬ 
servational effects. We first convert the linewidth from the fit to a velocity dispersion, 
cr obs = Au(8ln2) -1 / 2 . The data are discretely sampled along the spectral axis and thus 
limited by the spectral resolution of the VLA observations. The effective velocity resolution 
is <7 res ~ 0.6 km s -1 . Then 


Oiine ~ 



(2) 


is the “true” velocity dispersion along the line of sight. These values of the velocity dispersion 
averaged over the cprops clump boundaries are listed in Table [5] 


We further calculate the one-dimensional thermal linewidth calculated from the kinetic 
temperature, 

k-B Tr 


(Tj 1 — 


+NH 3 U^H 


( 3 ) 
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At; (km s 1 ) r o 


Fig. 13.— Plots of the results of NH 3 spectral line fitting. The circles are (Xr ed ) 1//2 - 
weighted averages over clumps defined by cprops, and stars are (y^ ed ) 1//2 -weighted aver¬ 
ages over whole IRDCs. Symbols are colored by 1RDC: G010.74-00.13, G022.56-00.20, 
G024.60+00.08, G028.23-00.19, G031.97+00.07, G032.70-00.30, G034.43+00.24, G035.39- 
00.33, and G038.95-00.47 are maroon, red, orange, yellow, bright green, dark green, cyan, 
blue, and magenta, respectively. The background black contours show the (x 2 ed ) 1//2 -weighted 
density of all lines of points in the sample (individual lines of sight and separated by distinct 
velocity components). Panels (a) and (b) also show red representative covariance ellipses 
from the spectral line fitting in the upper right. 


where /inh 3 = 17.03 is the mean molecular weight of ammonia and mn is the mass of a 
hydrogen atom. We can then calculate the uonthermal contribution to the velocity dispersion 











39 


via 



(4) 


We measure dispersions typically around 2 km s -1 , greater than the thermal linewidths, 
which are typically less than 0.1 km s -1 . The effect of this nonthermal dispersion is discussed 
in 


5.4 


5.3.2. Comparison to Previous Studies 


Taking the optical depth to be typically 5-10 for the majority of clumps, and the kinetic 
temperature to be 12-25 K for the majority of the clumps, we compare our results to similar 
studies. These values are in agreement with the dense clumps in G19.30+0.07 observed by 


Devine et al. 

(2011 

) with the VLA and the same spectral lines. Compared to the sample 

of 

Ragan et al. 

(20lTJ) , who also used the (1,1) and (2,2) lines observed by the GBT and 


the VLA, we find slightly higher kinetic temperatures in our study. They measured kinetic 
temperatures in the range of about 8 to 13 K, however their sample was selected to be devoid 
of star formation indicators and thus be in the earliest evolutionary phases. It is possible 
that our sample generally reflects a slightly later stage in IRDC evolution, in which star 
formation activity has increased and the gas is showing the affects of protostellar heating. 
Indeed, the clumps in our sample that are coincident with one or more 70 /jrn point sources in 
the Herschel images have kinetic temperatures 2-3 K higher on average than clumps without 
70 /jrn point sources. 


Our clumps are slightly different than dense cores and core candidates in Perseus ob¬ 
served by Rosolowsky et al. (2008) with the GBT, also with the same spectral lines. They 
found colder kinetic temperatures around 10-12 K and optical depths usually less than 5, 
but as high as 15 in the densest cores. Furthermore, they observed column densities around 
5 x 10 14 cm -2 with only the densest cores exceeding 10 15 cm -2 , while we observe column 
densities 10 14 -10 16 cm -2 and nearly all of our clumps have peak column densities of about 
3 x 10 15 cm -2 . We would expect higher column densities in our sample given that we selected 
IRDCs, which must have high enough dust column densities to be seen in extinction at MIR 
wavelengths. Their GBT observations of cores at a distance of 260 pc have comparable phys¬ 
ical resolution (0.04 pc) to our study (0.02-0.1 pc), so the differences cannot be explained by 
beam dilution effects. The higher column densities in our study can be partially explained 
by using a slightly different method of calculation that includes a correction for both parity 
states of the (1,1) (Friesen et al. 2009). This correction typically increases the column den¬ 
sities by less than a factor of two for the measured excitation temperatures. We therefore 
are observing slightly warmer gas and slightly higher column densities in IRDCs than lower 
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mass and apparently more quiescent star-forming environments like Perseus. This is not 
surprising if star formation activity is largely controlled by mass surface density; the average 
mass surface density in Perseus is approximately 90 M 0 pc -2 (Heiderman et al. 2010), but is 
approximately 150 M 0 pc -2 in IRDCs (based on masses and sizes reported by Simon et al. 
(2006b)). 


6. Discussion 

6.1. Kinematics and Previous Studies of Individual Sources 

All nine IRDCs show evidence of clumps and velocity substructure. The IRDCs are 
generally composed of many distinct clumps that sometimes can be grouped into distinct 
velocity components. These clumps also overlap spatially and show apparent interactions 
in position-position-velocity space, with star formation tracers (IR point sources, masers, H 
II regions, etc.) coincident with these overlapping sites. We also observe smoother velocity 
gradients along the most filamentary (high aspect ratios, A 0 > 3) clumps and across IRDCs 
as a whole. It is clear that the internal structure of varies from IRDC to IRDC, both in terms 
of the complexity of velocity substructure and in the prevalence of filamentary structures 
compared to the prevalence of globule structures. We discuss the characteristics of the 
individual IRDCs below. 


6.1.1. G010.74-00.13 


G010.74-00.13 shows a clear velocity gradient west to east across the IRDC. Channel 
maps of the NH 3 (1,1) emission are presented in Figure 14 At low velocity (vlsr ~ 28 km 


s _1 ), the NH 3 morphology appears filamentary, however it gradually transitions to a more 
globular morphology at higher velocity (ulsr ~ 30 km s~ x ). This structure appears to have 
two distinct velocity components with different morphologies. It is notable that both IR point 
(possible YSOs; see Figure [4]) are coincident with the filamentary component. There 


sources 


-i' 


m 


is an elevated velocity dispersion (peak in the second-moment maps of about 0.9 km s 
the center of the IRDC where these two velocity components overlap and may be colliding 
in this scenario. This location is also coincident with one of the protostcllar candidates. The 
majority of the NH 3 (2,2) emission is within the filamentary velocity component, whereas 
the CCS more closely follows the 30 km s -1 globular component. In G010.74-00.13, it may 
be that this cloud-cloud collision has triggered both the formation of the YSO and the high 
velocity dispersion, or it may be that the YSO is responsible for the high velocity dispersion 
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directly and its location is coincidental. 


6.1.2. G022.56- 00.20 


G022.56-00.20 is deconvolved into three clumps with cprops. The central velocity com¬ 
ponent (vlsr ~ 77 km s' 1 ) is located at the northeastern end of the IRDC, while the high 
(^lsr ~ 78 km s' 1 ) and low (wlsr ~ 75 km s' 1 ) velocity components have high aspect ratios 
(kL 0 ~3) and run nearly parallel to each other along the IRDC’s major axis. The 75 km s ' 1 
velocity component is located slightly to the north of the 78 km s ' 1 velocity component, 
though they overlap spatially, and the IRDC shows a high velocity dispersion of about 2 
km s ' 1 in the overlap region. It is noteworthy that the peak in integrated intensity, the 
peak in velocity dispersion, and an infrared point source are all coincident in the center of 
this overlap region, as can be seen in Figures [5] and 15 In Figure 16, the position-velocity 
diagrams taken across the IRDC parallel to the minor axis clearly show these major velocity 
components and their interaction. We identify two possible scenarios for G022.56-00.20: (1) 
this IRDC represents a collision of multiple velocity components, apparently triggering the 
formation of the YSO; or ( 2 ) the YSO is driving expansion of the molecular gas around it. 
The separation between the two primary components is < 5 km s ' 1 and they extend < 1 
pc, making either scenario plausible. 


6.1.3. G024.60+00.08 


G024.60+00.08 shows a clear gradient in the velocity field from about ulsr ~ 50.5 km 
s ' 1 in the west to about vlsr ~ 55 km s ' 1 in the east-southeast. Channel maps of the IRDC 
in Figure [17] do not clearly distinguish between the possibility of two distinct components or a 
single gradient across one major component. There is a velocity dispersion peak of about 1.2 
km s' 1 , cospatial with a large NH 3 clump and an IR point source near the center of the IRDC 
(see Figure [ 6 ]). Rathborne et ah (2007) identified two protostcllar condensations (bright, 
compact, millimeter cores with IR emission indicative of star formation) in G024.60+00.08 
with IR.AM Plateau de Bure 1.3 and 3 mm continuum and Spitzer images, one of which is 


the central IR point source. Cyganowski et ah (2008) identified an EGO in G024.60+00.08. 


The position of the EGO is marked in Figure [ 6 j and is also coincident with the protostellar 
candidate. 
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6.1.1 G028.23-00.19 


G028.23-00.19 is an apparent starless, dark cloud. The bright point source (Figure [T]) 
is an unrelated, foreground late-type star (Bowers & Knapp 1989). Rathborne et al.| ( |2006D 
made millimeter continuum maps of G028.23-00.19 with the 1RAM 30 m single dish telescope 
and observed three dense clumps at 11" resolution with masses ranging from 38 to 705 M 0 . 
The primary 1RDC has one filament, with at least 4 connecting filaments extending beyond 
our maps, but seen in 1R extinction. Sanhueza et al. (2013) found two distinct velocity 
components in CARMA observations of this IRDC, with 0.3 km s -1 , 1079 resolution. In our 
data, the position-velocity diagram shows that the velocity structure in NH 3 is a smoother 
gradient across the filament, from about ulsr ~ 79 km s -1 to ulsr ~ 82 km s -1 . It is unlikely 
that we are seeing a resolution effect, given that we have comparable velocity resolution and 
higher angular resolution than the CARMA study. The discrete clumps may be explained by 
a combination of chemical differentiation across the IRDC and/or missing the largest spatial 
scales in the CARMA data. There is also an ammonia velocity dispersion peak of about 0.9 
km s ^ 1 coincident with the optical depth peak. It is possible that the core in G028.83-00.19 
may be a deeply embedded YSO with an outflow exciting the SiO and CH 3 OH. However, 
this is unlikely given the lack of IR emission, the absence of a T K increase in the ammonia 


fitting results toward this source, and the detection of NH 2 D by Sanhueza et al. (2013). 
These characteristics are most consistent with a cold core. 


6.1.5. G031.97+00.07 


G031.97+00.07 has the most complex substructure in our sample. Rathborne et al. 


(2006) made millimeter continuum maps of G031.97+00.07 with the IRAM 30 m single dish 
telescope at 11" resolution and observed nine dense clumps with masses ranging from 151 to 


1890 M 0 . Wang et al. (2006) reported H 2 0 masers in G031.97+00.07 and Urquhart et al. 


(2009) identified an H II region in G031.97+00.07, all marked in Figure[ 8 ] The IR morphology 
has a long, thin filamentary structure leading to a higher contrast globule neighbored by at 
least 3 YSOs with H II regions and masers. The IRDC itself is part of a much larger 
molecular complex seen in 13 CO in the BU-GRS with IR dark filaments extending beyond 
our observations. 

Seen in Figure | 8 j the NH 3 emission closely matches the IR contrast, and it is deconvolved 


into 21 distinct clumps. As seen in the position-velocity slices in Figure 18, the structures 
are a mix of filaments and globules, and tend to fall in one of two distinct velocity ranges: 
Vlsr ~ 92-99 km s~ x and Vlsr ~ 97-102 km s -1 . The majority of the emission is in the 92-99 
km s ^ 1 velocity range, however the 97-102 km s -1 velocity range is coincident with at least 
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one YSO. The filamentary structures also show velocity gradients spanning about 2 km s _1 
along their major axis, as seen in Figure [19] This IRDC also has the strongest CCS emission 
in our sample, with a peak signal-to-noise ratio of about 22. G031.97+00.07 may be a region 
where several weaker filaments extending tens of parsecs are feeding molecular gas and star 
formation is progressing most quickly at the collision points. The morphology and velocity 


structure of this IRDC is consistent with the “hub-filament structure” described in Myers 


(2009) and Li et al. (2013), in which gas flows along the filaments to a central hub where it 


feeds star formation. 


Battersby et al. (2014b) recently studied G031.97+00.07 (called G32.02+0.06 in their 


work) in the NH 3 (1,1), (2,2), and (4,4) transitions with the VLA and discussed the environ¬ 
ment of the larger molecular complex. They observed two pointings, one toward the active 
region at the north end of the IRDC, and one quiescent region near the south end of the 
IRDC (their pointings partially overlap our combined maps, but do not cover the central 
portion of the IRDC). They found dense parsec scale filaments of 10-100 M 0 with dense 
cores less than 0.1 pc in size, in agreement with our findings. The authors reported that the 
dense cores were virially unstable to gravitational collapse, and that turbulence likely set 
the fragmentation length scale in the filaments. In the quiescent region, they observed the 
two distinct velocity components that continue along the IRDC toward the more complex 
active region in our data. Finally, they note the existence of at least three bubbles all seen 
in Spitzer , Herschel, and the BU-GRS that are likely older H II regions from previous gen¬ 
erations of massive stars, and may have compressed the molecular gas to form and/or shape 
the IRDC and trigger more recent massive star formation. 


6.1.6. G032.70- 00.30 

G032.70-00.30 has the weakest NH 3 emission in our sample. As shown in Figure [5] the 
(2,2) line is only marginally detected, and there is no detected CCS. Weak IR point sources 
indicate protostellar candidates at both ends of the IRDC. There is a velocity gradient from 
^lsr ~ 89 km s -1 to i+sr ~ 91 km s _1 across the filament from southeast to northwest, with 
the highest velocity dispersion of 0.8 km s -1 near the southwestern end. 


6.1.7. G03143+00.24 


Wang et al. 

(2006 

) reported H 2 0 masers in G034.43+00.24 and 

Urquhart et al. 

(2009) 


identified an H II region in G034.43+00.24, all marked in Figure [Id] Rathborne et al. 
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(2005) observed the millimeter/submillimeter continuum in G034.43+00.24 with IRAM, the 


James Clerk Maxwell Telescope (JCMT), and the Caltech Submillimeter Observatory (CSO). 
They identified three compact clumps of several hundred solar masses each. The SEDs 
stretching from the millimeter to the IR indicated high luminosities consistent with YSOs 
of ~10 Mq each. Moreover, Rathborne et al. (2005) also observed HCN, CS, and SiO in 
G034.43+00.24 with IRAM and CSO. The large line widths of ~10 km s -1 in HCN and 
CS along with the detection of SiO indicated outflows and shocked gas, further evidence of 
ongoing star formation. Cyganowski et al. (2008) identified two EGOs, marked in Figure 


10 Sanhueza et al. (2010) observed this region in multiple molecular gas tracers with the 


APEX 12 m telescope, the Nobeyama Radio Observatory (NRO) 45 m, and the Swedish- 
ESO 15 m Submillimeter Telescope (SEST). They found 4 molecular cores with velocity 
profiles indicative of outflows and large scale infall toward the most massive core, further 
strengthening the evidence for ongoing massive star formation in this IRDC. 


G034.43+00.24 has a complex substructure. The overall shape is very filamentary 


see 


Figure 10), and it is a portion of the Giant Molecular Filament GMF38.1-32.4a identified 


by Ragan et al. (2014). The IRDC is known to have three YSOs with masers, EGOs, and 


millimeter cores along the primary filament, with an additional YSO at the northern end of 
the IRDC. Shown in Figure |20j The velocity dispersion is elevated to over 2 km s -1 at two 
of the YSO positions. There is an overall gradient from the southwest to the northeast. The 
major filament has a velocity gradient from the west to the east, leading into a gradient from 
south (wlsr ~ 56 km s -1 ) to north (ulsr ~ 60 km s -1 ) along the weaker filament toward 
the northern YSO. This gradient may be the result of gas flowing along the larger GMF, or 
may be composed of multiple unresolved velocity components. 


6.1.8. G035.39-00.33 

G035.39-00.33 is a long, filamentary, high IR contrast IRDC with protostcllar candidates 
(see Figure [IT]). The velocity structure shows two distinct velocity components: a r^sR ~ 43 
km s" 1 velocity component at the northern end of our observations, and a i+sr ~ 45 km s^ 1 
velocity filament with a gradient from the southeast to the northwest. The two components 
are well separated in velocity space, and position-velocity slices do not show any apparent 
interaction between the two. We observe high velocity dispersion around most of the YSOs, 
only one of which is coincident with the lower velocity component. Our GBT data, which 
cover much more of the IRDC and have higher velocity resolution, show that the velocity 
gradient in the 45 km s^ 1 component continues along the northern extent of the filament. 
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Jimenez-Serra et al. (2010) observed G035.39-00.33 in the SiO J=2-l line (a shock 


tracer) with the IRAM 30 m, and found both bright, compact emission with broad linewidths 
and weaker, extended emission with narrow linewidths. The compact emission was consistent 
with protostellar outflows, and they proposed three explanations for the extended emission: 
(1) a collection of low mass outflows throughout the 1RDC, (2) a recently processed, more 
massive outflow whose energy was distributed throughout the 1RDC, or (3) the 1RDC was 
recently assembled from collisions of smaller clumps. They detected three filaments with 
distinct velocities at approximately 44.1, 45.3, and 46.6 km s _1 . Jimenez-Serra et al. (2014) 
further found velocity gradients of 0.4 to 0.8 km s -1 pc” 1 from north to south along the fila¬ 
ments, as seen in multiple transitions of 13 CO and C 18 0 with the IRAM 30 m and the JCMT. 
They further found that dense cores were preferentially located at the intersection of the (pos¬ 
sibly interacting) filaments. Henshaw et al. (2014) observed N 2 H + (1-0) in G035.39-00.33 
with the Plateau de Bure Interferometer (PdBI) array and noted that the local kinematics of 
the filaments were affected the dense cores, either through accretion or expanding envelopes. 


Hernandez et al. (2012) found that this parsec-scale region was consistent with being in 


virial equilibrium without requiring support from magnetic fields, and that the timescale for 
achieving equilibrium was less that the estimated 1.4 Myr for formation from two converging 
flows. 


While these studies focused primarily on the narrower portion of the IRDC north of our 
VLA maps, our larger GBT maps cover the region in these studies, and we also detect these 
three velocity components in NH 3 (though the 46.6 km s” 1 component is seen only weakly). 
The two velocity components we observed in our combined maps are the extensions of two of 
these filaments to the south. Most of the emission we detected here is part of the 45 km s” 1 
velocity component, and this southern region is both wider and darker (in MIR extinction) 
than the northern portion of the IRDC. We suggest that this IRDC may also represent an 
example of the hub-filament structure, in which gas is flowing along the parsec-scale 45 km 
s” 1 filament to the region we studied, where it feeds the ongoing star formation traced by 
IR point sources. 


6.1.9. G038.95-00.47 

G038.95-00.47 is a globular IRDC with long, low IR contrast filaments extending from a 
few parsecs away from it to the east. Seen in Figure [12} it also contains at least 2 protostellar 
candidates along the western edge, and neighbors two IR bubbles to the southeast and 
southwest. The velocity field shows a clear gradient from the southwest to the east (ulsr ~ 
40.5 km s” 1 ) and north (ulsr ~ 43.5 km s _1 ). This IRDC is also consistent with the “hub- 
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filament structure” described in Myers (2009) and Li et al. (2013). 


Xu et al. (2013) studied G038.95-00.47 in relationship to its environment, specifically 


the two adjacent H II regions G38.91-0.44 (N74) and G39.30-l.04 (N75). They used 1.4 
GHz radio continuum and single dish 12 CO, 13 CO, and C 18 0 J=l-0 observations to probe 
the physical characteristics of the bubbles and used GLIMPSE data to identify YSOs in the 
region. They found a molecular clump associated with G038.95-00.47 that showed a sharp 
integrated intensity gradient toward the H II regions consistent with it being compressed by 
their expansion. They further identify an overabundance of YSOs in the IRDC and propose 
that star formation was triggered via radiatively driven implosion. 


6.2. Mass Estimates 

We estimate the H 2 masses for IRDCs from three independent methods: dust extinction, 
dust thermal emission, and 13 CO emission, described below. Total IRDC masses for our 
sample for all three methods are shown in Figure |2TJ with a comparison discussed in 
All masses were measured by integrating over the footprint of our combined GBT+VLA 
observations. 


16.2.4 


6.2.1. IR Extinction Mass Estimates 


For the dust extinction method, we continue to follow the analysis presented by Butler 

The mass surface 


& Tan (2009) used to calculate IR extinction maps described in 14.2 


density is determined from the optical depth by 

y _ D 

— —• 
K A 


(5) 


where is the dust opacity per gas mass at wavelength A. Our choice of dust opacity from 
the literature will affect the mass estimate by as much as factors of a few. For MIPS 24 /irn, 
we adopt Hu /,m = 13.3 cm 2 g ~ 3 from the dust model of Ossenkopf & Henning (1994) for 10 5 
year-old, 10 6 cm -3 dust with thin ice mantles and assuming a dust-to-gas mass ratio of 1 : 100 . 
This value has an uncertainty of about a factor of three given the full range of possible 
values, depending on the age, thickness of ice mantles, and number density of dust grains, 
which results in a factor of about three uncertainty in the mass estimates. The angular 
resolution of the 24 /irn images, and thus our extinction and mass surface density maps, is 
6 ", compared to the 3'.'2-5 "A resolution of the NH 3 maps. This is the closest resolution match 
of the three mass estimates, however it lacks any velocity information to disentangle multiple 
clumps along a single line of sight. 
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For the total mass, we simply sum over the IRDC: 


M m = D 2 n ^ 

IRDC 



IRDC masses estimated with this method range from about 280 M 0 to 3900 M 0 , excluding 
G034.43+00.24. In this IRDC, Mjr is a very weak lower limit on the mass because there is 
significant IR emission from IR point sources that prevents us from calculating the extinction 
across most of the IRDC. 


6.2.2. BU-GRS 13 CO J=1-0 Mass Estimates 


Next, we estimate the gas mass from the 13 CO J=l-0 BU-GRS data cubes. Because the 
spectral line data have velocity information, we restrict measurements from this emission to 
velocity ranges with NH 3 (1,1) emission greater than 5a. While the velocity resolution of 0.2 
km s _1 is in fact better than our 0.6 km s ” 1 resolution NH 3 data cubes, the beam FWHM is 
46", so we are not able to resolve individual clumps. G010.74-00.13 is outside the coverage 
of the BU-GRS, and so is not included in this analysis. 

To calculate the 13 CO J=l-0 column density, we assume that the gas is in LTE, fills 
the beam, and is optically thin, and thus use the standard equation: 


1V( 13 C0) 


3h 


U3CO 


8x 3 U\ 3GO S ( J u ) // ] 3 (jq SuSK ^nuclear 
E u , 13CO 


x exp 


fcBD ex ,13CO 


- B,13CO 


dv, 


(7) 


where z+ 3 co = 110.201353 GHz, S(J U ) = J U /(2J U + 1) is the statistical weight of the upper 
level, J u = 1 for the upper level, Hi 3 co = 0.112 debye is the electric dipole moment, g u = 
2 J u + 1 is the rotational degeneracy, g k = 1 is the K degeneracy for a linear molecule, 
^nuclear = 1 is the nuclear spin degeneracy, and E u ^ 3CO /kB = 5.29 K for the energy of the 
upper level. We assume a partition function of the form Zi 3 co ~ 0.38(T eX) i3co/K + 0.88) 
(Wilson et al. 2009) and an excitation temperature T eX! i 3 co = 10 K. Using a higher excitation 
temperature of 20 K would raise the mass estimates by approximately 40% compared to 10 K. 
We adopt a conversion factor iV(H 2 )/iV( 13 CO) = 5 x 10 5 from Simon et ah (2001) (assuming 


R( 12 CO/ 13 CO) = 45 ( Langer &; Penzias|l990 ) and X( 12 CO) = 8 x 10 5 ( Blake et~ah|l987 )), 
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and thus the column density of the total molecular gas is 


iVi3Co(H 2 ) = 4.2 x 10 20 


// ?B,13CO dv\ 

V K km s _1 / 


cm 


To calculate the total mass of the gas, we sum over the 1RDC and use 


Misco — liTnuD 2 Q -^1300 (H 2 ) 

IRDC 

D \ 2 f n \ / hsco 
kpc ) \ arcrnin 2 ) \K km s" 1 



( 8 ) 


(9) 


where n is the mean molecular weight in multiples of the proton mass (assumed to be 2.33), 
D is the distance to the region, hi is the solid angle occupied by the gas, and / 13 C 0 is the 
sum of Xb,i 3 CO over the IRDC in position and velocity. 


IRDC masses estimated with this method range from about 2090 M 0 to 36,000 M 0 . 
We note that our mass estimates from this method are systematically lower than those 


reported by Simon et al. (2006b) by factors of about two to six. However, the Simon 


et al. (2006b) analysis uses ellipses larger than the IRDCs themselves, sometimes including 


complex networks of filaments extending beyond our maps, to get total masses, and so they 
expect their own values to be overestimates. Both our mass estimates with this method and 
those of Simon et al. (2006b) are systematically higher than those derived from 1R extinction. 


We discuss this further in f 6.2.4 


6.2.3. BGPS Dust Emission Mass Estimates 



iVmm(H 2 ) = 


Si. 


12mm 


^mb [MTl H & 1.12mm B\ (Td ) 


( 10 ) 


where S \,\ 2mm is the flux density of the Bolocam emission, Q m b is the solid angle of the main 
beam (33" FWHM), // is the mean molecular weight, mn is the mass of a hydrogen atom, 
r i.i 2 mm is the dust emissivity per unit gas mass, and B\(T d ) is the Planck function for a 
dust temperature T d . We adopt «i.i 2m m = 0.012 cm 2 g -1 , again assuming a dust-to-gas mass 


ratio of 1:100 and using the dust model of Ossenkopf & Henning (1994) for 10 5 year-old, 10 6 
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cm 3 dust with thin ice mantles. The total mass summed over the IRDC is then 


= 41 i ^ 


M mm = fJtm H D 2 n ^2 

IRDC 

D \ 2 ( n 

Jy ) \kpcJ \arcmin 2 

12.85 Kl \ „ r 

1 M 0 . 


x exp 


T, 


( 11 ) 


Our choice of dust temperature is informed by the findings of Battersby et al. (2014a). They 


compared the dust temperature inferred from Herschel SEDs to the gas (kinetic) temperature 
measured from the NH 3 (1,1), (2,2), and (4,4) transitions with the VLA in G031.97+00.07 
(called G32.02+0.05 in their work). The dust temperatures, 11.6+0.2 K, were systematically 
lower than the gas temperatures, 15.2+1.5 K, however the two temperatures were generally 
correlated regardless of local star formation activity. We therefore adopt the median of 
the clump-averaged kinetic temperatures in Table [5] for each IRDC multiplied by the ratio 
11.6/15.2 as the dust temperature for each IRDC. We do not use our T K maps directly 
because of the resolution mismatch and the lack of kinematic information in the Bolocam 
data to separate multiple velocity components. IRDC masses estimated with this method 
range from about 740 M Q to 8900 M 0 . It should be noted that the resolution of the BGPS 
is poorer than that of the ammonia maps by a factor of a few, and also does not contain 
velocity information. 


6.2.4. Comparison of Mass Estimates 


The IRDC mass estimates from the three methods are shown in Figure 21 For the 


majority of our sample, the masses derived from dust emission are typically consistent with 
those estimated from the IR extinction within a factor of a few, well within our estimated 
uncertainties. We note that in G028.23-00.19 and G032.70-00.30 , M IR and M 13CO agree 
well, but M mm is well below the other estimates. These are both quiescent clouds, so it is 
possible that the typical dust temperature in these IRDCs is lower than the median clump- 
averaged kinetic temperature we used, and so the dust emission mass estimate is low. I 11 
G034.43+00.24, Mir is a very weak lower limit on the mass because of IR emission from IR 
point sources, as discussed above. 


In every case, M^co is greater than M mm by at least factors of a few, and greater than 
an order of magnitude in G024.60+00.08. However, these two estimates are clearly well 
correlated. The 13 CO <7=1-0 emission is probably optically thick, and may miss gas mass 
because of CO depletion, but also includes the contribution of the lower density molecular 
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envelopes surrounding the IRDCs. Hernandez & Tan (2011) compared the mass surface 


density estimates from 8 /irn extinction to that from the BU-GRS, taking care to subtract 
the contribution from the molecular envelope surrounding the IRDC itself (their study was 
performed in G034.43+00.24 and G035.39-00.33, though selected slightly different portions 
of the IRDCs that we studied here). Their mass estimates were also well-correlated though 
systematically offset, and they estimated the mass surface density of the envelopes to be 
about 50% higher that within the IRDC itself. Considering that the envelopes have a greater 
physical extent than the IRDCs, this effect can easily explain the difference in masses derived 
from dust and from 13 CO emission. 


6.3. Ammonia Abundance 


Total masses of the individual clumps are necessary to investigate their stability against 
gravitational collapse to form stars. We cannot use any of the three methods used above 
for the total IRDC masses because none of the datasets have both sufficient angular and 
velocity resolution to be used on individual clumps. It is straightforward to calculate a total 
NH 3 column from the spectral line fit, as shown in Appendix [Bj We then must assume a 
fractional abundance, A(NH 3 ), to convert from NH 3 column density to H 2 column density. 
Our ability to compute the total mass of molecular gas in individual clumps, and thus their 
gravitational stability, is limited by our knowledge of the ammonia abundance relative to 
molecule hydrogen. Previous studies of the fractional abundance of NH 3 in physically similar 
environments yield a range of values from 10~ 9 to ICR 7 (see for example Ragan et al. 2011 


Chira et al. 2013; Battersby et al. 2014a), and may depend on galaetocentric radius, IRDC 


age, and physical environment (Battersby et al. 2014a). The uncertainty in this conversion 
factor will translate linearly to the uncertainty in clump masses, and so will affect our analysis 
of gravitational stability. 

Because of the different resolutions and the fact that our ammonia column and mass 
measurements are restricted to the dense gas, we cannot probe how the abundance varies 


with environment. However, as shown in Figure 21, adopting a value of X(NH 3 ) = ICR 8 
for the ammonia abundance relative to hydrogen typically gives good agreement with dust 
emission and IR extinction estimates within a factor of about three. We therefore assume a 
nominal value of ICR 8 for calculating the masses of clumps, M c \, however the uncertainty of 
any individual clump mass is likely at least a factor of three. 
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6.4. Gravitational Stability 


When calculating the parameters concerning gravitational stability of clumps, we im¬ 
plicitly assume that all filamentary structures are elongated completely within the plane of 
the sky. We discuss the effects of three-dimensional geometry in Appendix [C] assuming a 
random distribution of orientations in space. We cannot correct values for individual clumps 
with existing observations, but the typical corrections for these values are less than a factor of 
2 and often close to 1. The projection effects will be quite large for clumps elongated nearly 
along the line of sight, but we expect that only a small fraction of clumps have this property. 
Values reported in this paper are presented without correction for projection effects. 


From the mass and size of a clump, we can calculate the mean mass density po, and 
thus the spherical free-fall time, 


^ff,sph 

-^eff 


3n 

32Gp 0 
3/2 'M, 



= 16.6 


or the filament free-fall time for a cylinder collapsing along its axis, 


( 12 ) 


^ff,cyl \j gdlg^fFjSph 


= 13.5An 


RcS \ 
pc ) 


3/2 


M r ] 


M, 


0 


- 1/2 


Myr, 


(13) 


where A 0 is the initial aspect ratio of the filament (Pon et al. 2011). We adopt the current 


aspect ratio of the clumps for Aq, implicitly assuming that these structures are either young 
or have not evolved substantially since their formation, so our values of tff, cy i will be lower 
limits. Spherical free-fall times for these clumps range from 2 x 10 4 to 2 x 10 5 years, and 
cylindrical free-fall times range from 2 x 10 4 to 10 6 years. These timescales are consistent 
either with the expectation that these IRDCs are still in the earliest phases of star formation 
or that they are supported against free-fall collapse. 

The first models of filamentary structures collapsing along characteristic length scales 


was presented by Chandrasekhar & Fermi (1953). They determined that a cylindrical cloud 


of incompressible fluid with an axial magnetic field separated into collapsing fragments whose 
length was initially 12 times the diameter of the cylinder. Following collapse and redistri¬ 
bution of mass, the length to diameter ratio shrinks to 4.5 if the average density remains 


unchanged, or 2 if there is no magnetic field. Schneider & Elmegreen (1979) identified fila¬ 


ments in optically dark nebulae and measured ratios of core separation to filament diameter 
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of 3 ± 1, consistent with these predictions. The aspect ratios of approximately two-thirds 
of the clumps in our sample are between 1 and 2 , but the most filamentary of our clumps 
have aspect ratios greater than 3, consistent with these predictions and studies. 


It is common to calculate the virial mass and virial parameter for spherical clumps via 


®vir,sph 


A/ V .jr_sph 

Md 


GM cl 


(14) 


using the total clump mass determined from NH 3 , M c \, the line of sight velocity dispersion 
u line, and the effective radius. A virial parameter less than 1 indicates a spherical clump prone 
to collapse, and greater than 1 is resistant to collapse. This formulation of the parameter 
assumes a uniform density sphere, which is unlikely to be an accurate description. If the 
clump is spherical but the mass density has a power law distribution, i.e. p oc r~ p , then 
the gravitational potential energy, and thus the virial parameter, will change by a factor of 
(15 — 6p)/(15 — 5 p). For example, the virial parameter will be a factor of 9/10 as large if 
p = 1 or 3/5 as large if p — 2. Power law density distributions can then reduce the virial 
parameters by about a factor of 2 , and so clumps will be more prone to collapse. 


To calculate a virial parameter for filamentary clumps, we use the equilibrium relation 
from Chandrasekhar & Fermi (1953) for an infinite, isothermal cylinder: 


2 ( 7 - 1 ) 



-G 


( A4 v i r , C yl \ 

V ^“) 


o. 


(15) 


We take 7 = 7/5 for the adiabatic index of diatomic gas, (14/L) is the internal kinetic energy 
of the clump per unit length, ( B/L ) is the magnetic energy per unit length, and (M vir)Cy i/L) 
is the cylindrical virial mass per unit length. We ignore the contribution from magnetic fields 
for now, and assume the internal kinetic energy is turbulence dominated since the linewidths 
are largely nonthermal, such that 

L4 _ 3 f Af v i r C yi^ 2 
L ~ 2 J ^ 



For cylinders of finite length, we note that the L & 2 R e g^G% by approximating the 
filament as an ellipse (as in cprops to calculate the aspect ratio) and noting that it has a 
projected area equal to 7rR^ s . We then recover an equivalent virial parameter for cylindrical 
clumps: 


12 

~5~ 



GM C \ 


a 


vir,cyl 


^^vir.cyl 

Md 


12 

25 


\J AqO! 


vir,sph* 


(17) 











53 


For the most filamentary clumps in our sample A 0 > 3 (with one being approximately 6 ), 
so a vil . iCy i ~ a v ir,s P h- For clumps with smaller aspect ratios a v ir,cyi < «vir,s P h, however this 
estimate is less reliable because the initial assumption of a cylinder of infinite length is 
certainly violated in this case. The median value of a v ir,s P h is about 1.5, and is about 1 for 

^vir,cyl • 

Given the factor of three uncertainty in the individual clump masses and the assump¬ 
tion of relatively simple geometry, we cannot confidently say which individual clumps are 
gravitationally bound; however, it is significant that across our large sample most clumps 
are approximately in equipartition of gravitational and kinetic energies. For clumps that are 
interacting and not in isolation, then the simple assumptions of virial equilibrium are not 
met and it becomes more difficult to determine if these structures are bound. However, if the 
clumps are completely collapsing, then the kinetic energy should be twice that of the grav¬ 
itational energy. Since the sample average kinetic and potential energies are approximately 
equal it is more likely that these clumps are prevented from collapsing by their internal 
motions, rather than that they are already collapsing. This is consistent with the fact that 
the star formation is not predominantly in the center of the clumps, but at the interfaces 
between clumps. 


We further investigate the role of turbulent support in these clumps by calculating the 
sound speed of the gas as given by Chandrasekhar & Fermi (1953), 


and then the Mach number, 


c, = 


M = 


lk B T K 

(18) 

pm H 

°NT 

(19) 


We take 7 = 7/5 as above, and list the sound speeds and Mach numbers for the clumps in 
Table [ 6 ] The Mach numbers are typically 3 to 5 and all greater than 1 (with the exception of 
one clump in G034.43+00.24 for which the linewidth is approximately the thermal linewidth). 
This implies that the clumps are dominated by nonthermal motions. These nonthermal 
motions may be dominated by turbulence that provides support against gravitational collapse 
of the clumps. It is also possible that high nonthermal motions arise from YSOs, for example 
from infall motions, which do not provide support. However, it is unlikely that gas motions 
directly associated with YSOs can explain the high nonthermal motions across the entirety 
of the IRDCs, so turbulence is likely the dominate cause and thus an important source of 
support for these clumps. 
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Tablc 6. Clump Stability 


M cl 

(Mg) 

tff,sph 

(Myr) 

^ff,cyl 

(Myr) 

-^vir,sph 

(Mg) 

^^virjCyl 

(Mg) 

^vir,sph 

^vir,cyl 

c s 

(km s' 1 ) 

M 

B CI 

(mG) 

409 

0.06 

0.10 

217 

152 

0.53 

0.37 

0.26 

3.9 

1.37 

454 

0.05 

0.05 

171 

91 

0.38 

0.20 

0.27 

3.5 

1.88 

112 

0.06 

0.10 

59 

41 

0.53 

0.36 

0.26 

2.5 

0.87 

348 

0.06 

0.12 

195 

149 

0.56 

0.43 

0.27 

3.8 

1.36 

486 

0.05 

0.11 

287 

238 

0.59 

0.49 

0.29 

4.4 

2.09 

162 

0.04 

0.05 

264 

158 

1.63 

0.97 

0.28 

5.6 

1.95 

103 

0.04 

0.10 

257 

207 

2.50 

2.01 

0.28 

5.5 

1.33 

174 

0.05 

0.06 

271 

155 

1.56 

0.89 

0.31 

4.5 

1.28 

243 

0.07 

0.09 

136 

80 

0.56 

0.33 

0.26 

3.2 

0.88 

63 

0.06 

0.06 

40 

21 

0.63 

0.33 

0.27 

2.2 

0.72 

88 

0.05 

0.05 

119 

63 

1.36 

0.72 

0.30 

3.5 

0.99 

6 

0.04 

0.10 

3 

3 

0.52 

0.46 

0.24 

1.2 

0.67 

14 

0.05 

0.06 

10 

6 

0.70 

0.42 

0.28 

1.4 

0.59 

199 

0.03 

0.04 

156 

91 

0.78 

0.46 

0.26 

4.5 

2.34 

80 

0.04 

0.04 

85 

46 

1.06 

0.57 

0.26 

3.7 

1.30 

99 

0.04 

0.05 

199 

128 

2.01 

1.30 

0.25 

6.0 

1.73 

21 

0.03 

0.05 

41 

26 

1.93 

1.24 

0.26 

3.3 

1.05 

24 

0.05 

0.11 

82 

67 

3.45 

2.81 

0.27 

4.1 

0.72 

41 

0.03 

0.12 

74 

74 

1.81 

1.81 

0.25 

4.2 

1.31 

42 

0.03 

0.04 

32 

20 

0.77 

0.48 

0.26 

2.9 

1.84 

28 

0.02 

0.03 

58 

38 

2.09 

1.36 

0.24 

5.0 

2.30 

871 

0.04 

0.05 

212 

118 

0.24 

0.13 

0.28 

3.7 

2.89 

103 

0.05 

0.08 

36 

24 

0.35 

0.24 

0.25 

2.2 

1.18 

116 

0.05 

0.05 

631 

336 

5.43 

2.89 

0.31 

7.5 

1.14 

454 

0.05 

0.07 

301 

179 

0.66 

0.39 

0.30 

4.1 

1.62 

408 

0.06 

0.09 

319 

210 

0.78 

0.51 

0.30 

4.2 

1.45 

59 

0.04 

0.05 

39 

24 

0.65 

0.41 

0.29 

2.5 

1.37 

59 

0.04 

0.10 

87 

71 

1.46 

1.20 

0.29 

3.5 

1.16 

551 

0.05 

0.06 

568 

335 

1.03 

0.61 

0.29 

5.8 

1.95 

104 

0.04 

0.04 

225 

129 

2.17 

1.25 

0.29 

5.5 

1.77 
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Fig. 14.— Channel maps of the main hyperfine component of the NH 3 (1,1) emission in 
G010.74-00.13, shown in a logarithmic scale from 5 to 70 mJy beam -1 . Black contours show 
the integrated NH 3 (1,1) intensity and overlapping GBT and VLA footprint, same as the 
white contours in Figure [|J Channel velocities, vlsr, are labeled in the upper right corner 
of each panel in km s -1 . The 70 pm infrared point sources are marked by stars in panel (a) 
(red in the online figure). There is a clear gradient from the filamentary structure running 
north to south into the globular structure extending to the southeast. The largest linewidth 
is at the center of the IRDC, coincident with both structures and with an embedded infrared 
point source. 
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G022.56-00.20 GBT+VLA NH 3 (1,1) 

Velocity Field (km/s) Velocity Dispersion (km/s) 
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Right Ascension (J2000) 


Fig. 15.— Moment 1 and moment 2 maps showing the velocity field ((no p)) and velocity 
dispersion ((v^ ^)), respectively, from the NH 3 (1,1) emission in G022.56-00.20. Clump de- 
convolution with cprops identifies three clumps: one filament running along the northwestern 
edge, one globule at the northeastern end, and one filament running along the southeastern 
edge. The gradient in the velocity velocity field is explained by the velocity offset between 
the two filaments, and the elevated velocity dispersion through the center of the IRDC is 
coincident with the overlap of the two filaments. The peak in the velocity dispersion, the 
peak in the integrated NH 3 emission, and an infrared point source are all coincident at the 
center of the IRDC. 
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G022.56-00.20 
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Fig. 16.— Position-velocity diagrams in (a) through (d), with white contours at 2a, 4a, and 
8a. Slices for the diagrams are labeled and overplotted on the integrated intensity map of 
G022.56-00.20 in (e) in grayscale with a linear stretch. The black contours for the integrated 
intensity and combined GBT and VLA footprint are the same as the white contours in Figure 
[5| The angular offsets in the diagrams increase from east to west. The three most distinct 
structures in each diagram are the main and two inner satellite hyperfine components of 
the NH 3 (1,1) kne. The velocity offset between the two primary filaments is most apparent 
through the center of the IRDC. 
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Right Ascension (J2000) Right Ascension (J2000) Right Ascension (J2000) Right Ascension (J2000) 


Fig. 17.— Channel maps of the NH 3 (1,1) emission as in Figure 14 but for G024.60+00.08, 
with contours and held of view matching Figure [6j and only one 70 /ini source. There is 
a clear gradient from the northwest to southeast along the filamentary structures. These 
filaments extend further in the infrared extinction maps beyond the coverage of the NH 3 
observations. Deconvolution with cprops identifies multiple clumps with offset velocities 
along these filaments. The largest linewidth is at the center of the IRDC, coincident with 
both structures and with an embedded infrared point source. 
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G031.97+00.07 
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Fig. 18.— Position-velocity diagrams in (a) through (f) as in Figure 16, but for 
G031.97+00.07 with white contours at 2 a, 4cr, 8a, and 16er. Slices and the integrated 
intensity map are shown in (g) also as in Figure 16, but with contours taken from Figure [8] 
G031.97+00.07 contains a complex substructure of globules and filaments at offset velocities 
and with internal velocity gradients. 
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G031.97+00.07 GBT+VLA NH,(1,1) 
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Fig. 19.— Moment 1 and moment 2 maps showing the velocity field ((no i))) and velocity 
dispersion respectively, from the NH 3 (1,1) emission in G031.97+00.07. Decon¬ 

volution with cprops reveals a complex substructure with 21 globules and filaments. The 
majority of the emission is in a lower velocity component around 92-99 km s -1 , though a few 
clumps form a higher velocity component around 97-102 km s -1 . These velocity components 
are cospatial, as shown by the elevated values in the velocity dispersion map. There are 
additionally velocity gradients along the substructures with higher aspect ratios, such as the 
filament in the middle of the IRDC. 
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G034.43+00.24 GBT+VLA NH,(1,1) 


Velocity Field (km/s) Velocity Dispersion (km/s) 
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Fig. 20.— Moment 1 and moment 2 maps showing the velocity held ((foi))) and velocity 
dispersion ((i?^ ^)), respectively, from the NH 3 (1,1) emission in G034.43+00.24. The IRDC 
as a whole is very filamentary, but deconvolution with cprops reveals clumps ranging from 
globule to filamentary within the larger filamentary structure. The southern half of the 
IRDC contains three strong protostellar candidates. The two with known EGOs, indicative 
of molecular outflows, are coincident with the peaks in the velocity dispersion. An overall 
velocity gradient in the IRDC is seen from the western edge to the eastern edge in this part 
of the IRDC, then continues along the fainter filament to the YSO in the north. 
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Mass (Af 0 ) 


Fig. 21.— Total molecular gas masses for our IRDC sample using three different method of 
estimation, along with the mass inferred from NH 3 assuming a relative abundance of 1CF 8 
compared to molecular hydrogen. The error bars represent a factor of three variation in 
this abundance. Values of Mir are the masses estimated from dust extinction, i.e. 24 /mi 
contrast. Because of the infrared emission from protostcllar candidates, IR bubbles, etc., 
the dust extinction method is not viable everywhere. The mass estimates may therefore be 
lower limits, particularly for G034.43+00.24. Values of M mm are the masses estimated from 
dust emission at 1.12 mm, observed in the BGPS. Values of Mi 3 co are masses estimated 
from 13 CO J=l-0 in the BU-GRS, except for G010.74-00.13, which was not covered in that 
survey. The emission is probably optically thick, but also includes the contribution of the 
lower density molecular envelopes. The masses estimated from the BGPS and the BU-GRS 
data are typically well-correlated, but the BU-GRS masses are higher by factors of a few. 












63 


Table 6—Continued 


M cl 

(Me) 

^ff ,sph 

(Myr) 

^ff,cyl 

(Myr) 

^^virjSph 

(M s ) 

^^virjCyl 

(M s ) 

^vir,sph 

^vir,cyl 

Cs 

(km s' 1 ) 

M 

B C r 

(mG) 

93 

0.04 

0.08 

300 

219 

3.21 

2.35 

0.31 

5.7 

1.40 

96 

0.03 

0.07 

91 

80 

0.94 

0.83 

0.25 

4.5 

2.63 

26 

0.10 

0.11 

383 

212 

14.81 

8.20 

0.29 

6.3 

0.28 

715 

0.04 

0.06 

895 

560 

1.25 

0.78 

0.30 

7.2 

2.71 

109 

0.04 

0.11 

95 

84 

0.87 

0.77 

0.26 

3.7 

1.54 

2759 

0.05 

0.04 

1082 

536 

0.39 

0.19 

0.31 

5.9 

3.47 

654 

0.05 

0.08 

177 

119 

0.27 

0.18 

0.26 

3.6 

2.00 

2110 

0.05 

0.08 

817 

566 

0.39 

0.27 

0.29 

5.8 

3.39 

102 

0.03 

0.13 

209 

218 

2.04 

2.13 

0.25 

6.2 

1.78 

37 

0.04 

0.06 

104 

72 

2.84 

1.96 

0.30 

4.2 

1.14 

826 

0.04 

0.05 

786 

436 

0.95 

0.53 

0.30 

6.5 

2.69 

219 

0.05 

0.09 

120 

91 

0.55 

0.41 

0.30 

3.1 

1.59 

17 

0.03 

0.03 

9 

5 

0.54 

0.29 

0.23 

2.0 

1.32 

79 

0.04 

0.12 

54 

50 

0.69 

0.63 

0.27 

2.9 

1.31 

14 

0.05 

0.08 

9 

6 

0.65 

0.44 

0.25 

1.5 

0.54 

24 

0.10 

0.17 

222 

154 

9.27 

6.43 

0.31 

4.5 

0.27 

59 

0.04 

0.07 

216 

149 

3.65 

2.51 

0.31 

5.1 

1.06 

12 

0.06 

0.09 

25 

16 

2.12 

1.36 

0.25 

2.4 

0.39 

96 

0.04 

0.04 

670 

342 

6.97 

3.55 

0.35 

7.5 

1.41 

5 

0.06 

0.14 

6 

5 

1.35 

1.09 

0.24 

1.5 

0.31 

7 

0.04 

0.09 

71 

54 

10.43 

8.04 

0.26 

5.0 

0.53 

3 

0.04 

0.06 

4 

3 

1.50 

1.06 

0.26 

1.5 

0.48 

3 

0.04 

0.06 

43 

28 

16.55 

11.03 

0.26 

4.6 

0.42 

3 

0.04 

0.08 

5 

4 

1.98 

1.50 

0.24 

1.8 

0.44 

59 

0.04 

0.05 

449 

254 

7.56 

4.27 

0.37 

6.1 

1.06 

37 

0.07 

0.14 

103 

78 

2.79 

2.10 

0.26 

3.8 

0.48 

196 

0.06 

0.14 

489 

385 

2.49 

1.96 

0.32 

5.4 

1.01 

116 

0.06 

0.06 

305 

163 

2.63 

1.41 

0.31 

4.9 

0.97 

16 

0.05 

0.09 

59 

43 

3.65 

2.64 

0.27 

3.6 

0.58 

55 

0.04 

0.05 

156 

96 

2.86 

1.75 

0.33 

4.5 

1.41 
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In the absence of magnetic field strength measurements, we can estimate the critical 
magnetic held strength required to support uniform spherical clumps against gravitational 
collapse in the absence of turbulence, B cr , by taking the condition that 


3 GM\ 
5-Reff 


R 2 A 


effi 


( 20 ) 


so that 


B r 


1.02 x 10 


-4 


M r 


cl 


Mr. 


Ii 


eff 


-2 


rnG. 


( 21 ) 


If we postulate that the clouds are all supported entirely by magnetic pressure, the required 
magnetic held strengths to prevent gravitational collapse would be 13 /rG < B cr < 3.5 rnG, 


well matched to the typical magnetic held strengths measured in molecular clouds (Crutcher 


2012). These values of B cr are upper limits on the critical held strength since the support 


from turbulence and thermal motions has been ignored, so even smaller magnetic helds will be 
sufficient to stabilize the clumps. In fact, clumps with virial parameters greater than 1 require 
no magnetic held to stay supported against collapse. Typical magnetic helds strengths are 
therefore probably more than sufficient to support these clumps, though unnecessary for 
support given the supersonic turbulence. 


Comparing the free-fall times, the virial parameters, and the turbulent and magnetic 
support, the majority of clumps in these IRDCs appear to be supported against collapse. 
The observed nonthermal motions and typical magnetic held strengths in molecular clouds 
are sufficient to stabilize these clumps against self-gravity. It is likely that many of these 
clumps are then at least quasi-stable, and may be long-lived compared to the free-fall times. 
This may explain why we observe quiescent regions that are devoid of any signs of ongoing 
star formation, and why protostcllar candidates are typically found at the sites where clumps 
overlap. The clumps that are apparently interacting in position-position-velocity space are 
not evolving in isolation, and so collapse may be progressing in places where clumps are 
compressed by collisions with each other. 


6.J h l ■ Size-linewidth Relation 


We can place the clumps identihed in this study in the context of other molecular gas 
clumps in varying environments and physical characteristics by putting them on a plot of size 
versus linewidth or, equivalently, size and spherical virial mass. In Figure [22j red diamonds 
mark the dense clumps seen in 30 Doradus in the Small Magellanic Cloud in 13 CO J— 2-1 
with ALMA (Indebetouw et ah 2013). Green circles mark dense CS J— 2-1, 3-2, and 5-4 
clumps from a previous study of IRDCs by Gibson et al. (2009) with the 14 m FCRAO. 













Table 6—Continued 


M a 

(M 0 ) 

,sph 

(Myr) 

(Myr) 

-^vir,sph 

(Me) 

-^vir,cyl 

(M 0 ) 

^vir,sph 

^vir,cyl 

c s 

(km s' 1 ) 

M 

B cr 

(mG) 

116 

0.07 

0.16 

176 

144 

1.51 

1.24 

0.26 

4.1 

0.74 

76 

0.05 

0.06 

146 

83 

1.91 

1.09 

0.32 

3.7 

0.93 

3 

0.03 

0.05 

20 

13 

7.53 

4.87 

0.26 

3.5 

0.62 

58 

0.06 

0.10 

73 

50 

1.27 

0.87 

0.25 

3.3 

0.74 

249 

0.06 

0.07 

117 

71 

0.47 

0.28 

0.27 

3.1 

1.23 

12 

0.05 

0.07 

9 

6 

0.78 

0.51 

0.24 

1.7 

0.59 

55 

0.06 

0.06 

93 

52 

1.70 

0.95 

0.27 

3.6 

0.76 

6 

0.04 

0.07 

19 

13 

3.41 

2.33 

0.26 

2.7 

0.55 

126 

0.05 

0.05 

107 

57 

0.84 

0.45 

0.26 

3.6 

1.21 

106 

0.06 

0.08 

125 

79 

1.18 

0.75 

0.27 

3.7 

0.89 

187 

0.06 

0.07 

74 

42 

0.39 

0.23 

0.26 

2.6 

1.05 

115 

0.05 

0.05 

173 

94 

1.51 

0.82 

0.30 

4.2 

1.22 

1 

0.10 

0.16 

16 

11 

12.00 

8.22 

0.34 

1.8 

0.10 

26 

0.06 

0.10 

155 

104 

5.88 

3.93 

0.26 

5.3 

0.50 

71 

0.06 

0.09 

207 

129 

2.94 

1.82 

0.27 

4.8 

0.70 

7 

0.06 

0.05 

17 

9 

2.50 

1.27 

0.26 

2.2 

0.38 

9 

0.07 

0.09 

75 

45 

8.59 

5.18 

0.27 

4.1 

0.33 

55 

0.07 

0.11 

192 

130 

3.48 

2.36 

0.28 

4.7 

0.60 

61 

0.08 

0.07 

120 

61 

1.96 

0.99 

0.29 

3.3 

0.50 

68 

0.05 

0.05 

112 

64 

1.65 

0.94 

0.30 

3.6 

1.05 

6 

0.06 

0.07 

35 

19 

5.95 

3.29 

0.27 

3.0 

0.30 


Note. — Clumps appear in the same order as Table [5] 
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Caselli & Myers (1995) observed massive cores in Orion in 13 CO and C 18 0 J— 1-0, shown 


m 


blue crosses. Purple triangles mark galactic clouds listed by Heyer et al. (2009) as observed in 
12 CO and 13 CO 1-0 with the FCRAO, with their observed relation plotted in black (extended 
as a dotted line): M viriSp h = 300 M e R pc . Oka et ah (2001) observed clouds in the Central 
Molecular Zone in the galactic center in 12 CO 1-0 with the NRO 45 m telescope (plotted 
as orange X’s), and determined the relation M viTjSph = 2 x 10 4 M e R pc . Finally, Falgarone 


et al. (1992) observed quiescent molecular clouds in 13 CO 1-0 and 2-1 with the IRAM 30 m, 


shown as maroon squares. 


We additionally use the presence or absence of a 70 /im point source within the bound¬ 
aries of a clump to classify it as “protostcllar” or “starless,” respectively. This classification 
is listed in Table[5] Clumps from this work are shown as black stars in Figure[22] protostcllar 
clumps are filled, while starless clumps are empty. It is notable that the trend seen across 
different studies between radius and linewidth extends over more than 3 orders of magnitude 
in size, regardless of the molecular line tracers and observations used. We note that both 
classes of our clumps are along the same trend, though the protostellar clumps are typically 
larger in size and virial mass (i.e. linewidth). Recall that the thermal linewidths for even 
Tk = 50 K gas is only 0.12 km s _1 , and so the linewidths are largely nonthermal. 


Zhang et al. (2011) observed ammonia among other molecular tracers, and also saw 


that linewidths were elevated around protostellar candidates and were more elevated for 
more evolved sources. They further found that nonthermal velocity dispersions were su¬ 
personic. Urquhart et al. (2011) found that in massive cores the column densities, kinetic 


temperatures, and linewidths seen in ammonia are all correlated with each other and with 
bolometric luminosity, indicating that these values are driven largely by the central source. 
We also observe strong correlations between the column densities, kinetic temperatures, and 
linewidths in the clumps in our sample. 


Molecular clumps may be virialized but still have higher velocity dispersions (and thus 
virial masses) than Larson’s relation for galactic disk clumps if there exists some external 
pressure. Elevated linewidths may be a result of chaotic gravitational collapse, in which 
the turbulent motion arises from the collapse. An alternative explanation is pressure from 
the gravity of the larger IRDC as a whole. In this sample, we have already shown that 
the clumps have strong supersonic turbulence, which likely accounts for the relationship 
between the size and velocity dispersions (Ballesteros-Paredes et al. 2011). However, the 
velocity dispersions being elevated above the typical Milky Way relation requires the clumps 
to not be in equilibrium (as in the case of outflows), or have additional external pressure 
beyond that of the warm ISM. 
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Fig. 22.— A comparison of radius and spherical virial mass (a proxy for linewidth) for 
the clumps in our sample (using i? eff and M vir sph ) to other studies of molecular gas in 
different environments. Filled stars indicate clumps from our work that are determined to 
be protostellar due to the presence of a 70 /im point source within the clump boundary. Open 
symbols are starless clumps that lack a 70 /jrn point source. The black lines are relations 


from Oka et al. (2001) and Heyer et al. (2009) for their galactic center clouds and Milky 


Way disk clouds, respectively, extended to a greater range of cloud radii as dotted lines. See 
§6.4. 1| for descriptions of the different populations. 
















6.5. A(CCS)/A(NH 3 ) 



Suzuki et al. 

(1992 

) showed that a relative lack of CCS compared to NH 3 indicates 

a more evolved core. 

dirota et al. 

(2009 

) observed varying abundances of carbon-chain 


molecules relative to NH 3 in dark cloud cores and attributed the variation to the different 


evolution stages of the cores. Sakai et ah (2008) observed 55 millimeter cores in IRDCs with 


the NRO 45 m telescope and the Atacama Submillimeter Telescope Experiment (ASTE) 10 
m telescope. They observed N 2 H + , HC 3 N, CCS, NH 3 , and CH 3 OH at 18"-73" and 0.12 km 
s -1 to 0.5 km s' 1 resolution. They found that the [CCS]/[N 2 H + ] ratios of even 24 /mi dark 
cores were typically less than 1 and concluded that these cores were more chemically evolved 


than low-mass starless cores where the ratio is 2.6-3.2. Marka et al. (2012), however, found 


no variation in the column density ratio, A(CCS)/./V(NH 3 ), with evolutionary state of Bok 
globules, the low-mass dark cloud analogs to IRDCs. All of these studies targeted dense cores 


but were performed with single-dish observations. Devine et ah (2011) used the VLA and 


demonstrated that the spatial distributions of the NH 3 and CCS were highly anticorrelated 
in a single IRDC. Since these NH 3 and CCS transitions have similar excitation conditions 
(i.e critical densities), these studies show that neither molecule traces all of the gas. 


The spatial distributions of NH 3 and CCS in our sample are shown in Figures [4] 12 


Because of the low signal-to-noise in the CCS line, the GBT observations give a clearer 
impression the CCS distribution. The CCS generally is not detected over the full spatial 
extent of the NH 3 emission. In the scenario where the ratio of these abundances is a viable 
method for tracing evolutionary stage, we would expect the most evolved (i.e protostellar) 
cores would have a lower CCS abundance than starless cores. 


We investigate the possibility of chemical evolution of different clumps by comparing 


the column densities of NH 3 and CCS. Lai et al. (2003) calculate the column density of CCS 


in lower mass star-forming regions from the 22.34403 GHz line using an updated form of the 


relation from Suzuki et al. (1992): 


A(CCS) = 5.1 x 10 Tecs 


An, 


CCS 


X 


v km s' 1 
exp (E U} ccs/ [^B-^~ex,CCs]) 
exp (E U) ccs/ [&b -^ex, ccs] — 1) 


Zees 


cm 


-2 


( 22 ) 


where Auccs is the FWHM of the CCS line, Zees is the partition function (24 to 62 for low 
rotational transitions), E Ui ccs = 1-12 cm -1 is the energy level of the upper state, T eXi ccs is 
the excitation temperature of the line, and Tecs is the peak optical depth. The peak optical 
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depth is related to the peak main beam temperature by 


T, 


mb,CCS 


= $T, 


B,CCS 


$ \yf (^ex,CCS ) — J(T hg )] [1 — exp( —Tecs)] , 


(23) 


where T m b,ccs, ^b,ccs, and J(T) are the same as given in Appendix |B| but for CCS, Tb g = 
2.73 K is the background CMB temperature, and $ = 1 is the assumed beam-filling factor. 
Linewidths are difficult to measure for the CCS given that it is only detected with a signal- 
to-noise ratio greater than 5 along few lines of sight (refer to Figures |4p2|). We therefore 
take At’ccs from the ammonia velocity dispersion (second order moment) maps, U = 43, 
and T eX! ccs = 5 K. With these assumptions, we calculate column densities of approximately 
10 11 cm” 2 and as high as 10 12 cm -2 , similar to those observed by Suzuki et al. (1992) and 


Lai et al. (2003). Given the scatter in the points, the low signal-to-noise ratio of the CCS 


emission, and the assumptions made to obtain IV(CCS) we cannot determine the column 
densities along individual lines of sight to better than an order of magnitude. We therefore 
only use the clump averaged values of IV(CCS) to compare to the column density of the 

nh 3 . 

We expect from previous studies that CCS will be relatively depleted in the more chem¬ 
ically evolved gas. Without a direct probe of the age of the clumps, we use star formation 
activity as a proxy. Using the simple classification scheme mentioned in §6.4.1|, the presence 


or absence of a 70 //m point source within the boundaries of a clump are used to classify it as 
“protostcllar” or “starless,” respectively. We also consider the possibility that temperature 
and/or density may have an effect on the depletion of CCS onto dust grains. Figure 23 
shows a comparison of the ratio of the column densities of CCS and NH 3 against the kinetic 
temperature and the ammonia column density, averaged over clumps defined by cprops. 
The distinction between protostcllar and starless clumps in iV(CCS)/iV(NH 3 ) is not as clear 
as we would expect. Figure [23] shows marginal evidence of a trend for CCS to be relatively 
depleted compared to NH 3 in the highest column density gas, but one should be extremely 
cautious of interpreting this as a true correlation given the uncertainties involved. More 
sensitive high resolution studies of CCS, preferably including more than one spectral line, 
are necessary to say whether or not CCS is relatively depleted in certain environments in 
IRDCs. 


7. Conclusions 

In this study, we have mapped the NH 3 and CCS in nine IRDCs, combining single dish 
and interferometric data from the GBT and the VLA. From fitting the NH 3 spectral lines, we 
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Fig. 23.— A comparison of the ratio of the column densities of CCS and NH 3 against 
the kinetic temperature and the ammonia column density, averaged over clumps defined by 
cprops. Colors are as in Figure [13} Diamonds represent clumps that are coincident with a 
70 /jrn point source (protostellar), and circles represent clumps without such a point source 
(starless). The distinction between protostellar and starless clumps in At(CCS)/lV(NH 3 ) is 
not as clear as we would expect from previous studies. Given the scatter in the points, the 
low signal-to-noise ratio of the CCS emission, and the assumptions made to obtain IV(CCS), 
more sensitive high resolution studies of CCS, preferably including more than one spectral 
line, are necessary to say whether or not CCS is relatively depleted in certain environments 
in IRDCs. 



probed the physical conditions of these regions. We find typical values of physical parameters 
in agreement with other studies: kinetic temperatures around 12-25 K, linewidths around 
1-2 km s -1 , clump sizes from < 0.1 to 0.2 pc, and clump masses of tens to thousands of M 0 . 

The kinematics seen in NH 3 reveal a diverse set of configurations and substructure. 
It appears that no one description of internal structure works for all IRDCs. We observe 
IRDCs with gradients in the velocity field and velocity dispersion, IRDCs with discrete 
clumps possibly interacting with each other, and IRDCs that clearly have both of these 
scenarios occurring simultaneously in different locations. YSOs have a tendency to form 
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at the interface of different clumps. We find that these clumps are typically near virial 
equilibrium and can easily be supported by typical molecular cloud magnetic held strengths, 
so these structures may survive for multiple free-fall times. This result is consistent with 


turbulent core accretion models of massive and intermediate star formation (McLaughlin & 


Pudritz 1997 McKee & Tan 2003). The turbulent motions might also arise from the early 


phases of star formation and thus provide support against gravitational collapse only after 
some YSOs have been formed. 


At least three of these IRDCs are consistent with the picture of hub-filament structure, in 
which filaments feed molecular gas toward a cental, dense structure harboring star formation. 
All nine of these IRDCs have lower contrast IR extinction features that extend beyond our 
molecular maps, and so it is possible that all of these IRDCs are consistent with hub-filament 
structure. High resolution studies of other physical probes, including shock tracers like SiO, 
are necessary to further test whether the gas is in fact colliding in these regions. 

While this study has one of the largest samples of IRDCs studied at this level of detail, 
the variability of properties within this sample indicates that it is still not large enough to 
draw firm conclusions about the entire population of IRDCs. Aside from the dichotomy 
between globular and filamentary morphologies, there is also variation from more quiescent, 
slightly supersonic motion to highly supersonic motion. With evidence for at least two 
different classes of quite distinct morphologies and a possible hybrid classification but only 
10 sources we cannot determine with much significance which is most prevalent. A sample 
of a about 30 would show better if the categories are evenly populated and (about 15 in each 
category) or if one is more important and whether the hybrid classification is prevalent. 

The signal-to-noise ratio of the CCS line is typically low, and so deeper observations 
and observations of other CCS lines are necessary to better probe the possible chemical 
differentiation and evolution in IRDCs. 
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A. IR Extinction Maps 


First, a model of the IR background is generated using a 40' square median filter sampled 
every 24", and then smoothed by a 24" Gaussian filter. The median filter preferentially skews 
the model toward pixels lacking strong emission or extinction, provided the size of the filter 
is larger than the typical size of a dark cloud or other neighboring structures. Increasing 
the size of the filter also causes the model to be less local, so it is not beneficial to sample 
the model at the same pixel scale as the original images. The Gaussian smoothing removes 
unphysically discrete (sharp) features introduced by the discrete sampling of the model. 


The fraction of the IR emission that is foreground to the IRDC, /fore, is estimated assum¬ 
ing a galactic dust distribution, adopting a solar galactocentric radius of 8.4 kpc (consistent 
with our distance determination in and using the galactic longitude and distance for 


each IRDC. Following Butler & Tan (2009), we assume the hot dust follows an azimuthally 
symmetric distribution matching the surface density of Galactic OB associations, 


E ob oc exp 


R 


3.5 kpc 


(Al) 


where R is the galactocentric radius extending to 16 kpc. The values of / fore are listed in 
Table [5j and range from 0.058 to 0.328. If I\,o is the intensity of radiation just behind the 
IRDC and I\j is the intensity of radiation just in front of the IRDC, then they are related 
to the optical depth through the IRDC, T\, by 


h,i = 7 Ai0 exp(—r A ). (A2) 

Accounting for the foreground contribution, the intensity of our background model, Tyo.obs, 
is related to the true intensity just behind the IRDC by 


7a,0 (f /fore) /\,0,obs- (A3) 

Similarly, the observed intensity in the images, 7 Ajl)0bs , is related to the true intensity just in 
front of the IRDC by 


7a,1 7a, 1 ,obs /fore7A,0,obs- 


(A4) 
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Therefore, the optical depth through the IRDC is obtained by 


T\ = ~ 111 


7A,l,obs /fore7A,0,obs 

(1 /fore) l\,0,obs 


(AS) 


Battersby et al. (2010) make the point that extended sources in IRAC must be corrected 


for internal scattering in the instrument, and they include this correction for their analysis 


of the optical depth at 8 /irn. Engelbracht et al. (2007), however, find there is no need to 


apply such a correction for 24 /irn images from MIPS over a range of a factor of about five 
in background level, so we do not apply a scattering correction in our calculations. 


An optical depth map is then generated on the same pixel scale as the IR images. 
Point sources and other emission features appear as negative optical depth in the maps, 
and are thus excluded from analysis. Featureless portions of the IR image correspond to 
approximately 0 in the optical depth maps as expected. 


B. Ammonia Spectral Line Fitting and Physical Parameters 

The basic functional form of the fit to a single NH 3 inversion line is 

A T a ,M K ) = [J (Tex(J, K)) - J (T bg ( j, K))] [1 - exp (-r„(J, K ))], (Bl) 

where AT^ U (J, K) is the increase in corrected antenna temperature at the telescope from 
the line emission, rj m b is the main beam efficiency, $ is the beam-filling factor, T ex ( J, K) is 
the excitation temperature of the line, K) is the background temperature of the line, 

r u (J,K) is the frequency-dependent optical depth, and 


J(T) = + 
ks 


exp 


hv 

kwr 


-i 


(B2) 


is the Rayleigh-Jeans temperature for the Planck constant h, the Boltzmann constant k^, 
and frequency v (Ho & Town es|l98 3). We assume $ = 1 (i.e the brightness temperature, Tb, 
equals the main beam temperature, T m b). In principle the background temperature includes 
a contribution from the continuum temperature, Tc(J, K), such that Tb g (J, K) = Tc (./, K ) + 
T C mb, where T CM b = 2.73 K is the contribution of the cosmic microwave background. While 
we cannot exclude the possibility of weak continuum emission, we can set a limit on its 
contribution. The continuum emission is within the noise of 0 K over the majority of our 
observed area. The strongest possible continuum emission in our study is about 1 K, but is 
in a region in G034.43+00.24 with peak emission over 15 K and large linewidths such that 
the hyperfine components here are blended. This is then an upper limit on the continuum 
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emission taken from the lowest flux density channel between the peaks of the hyperhne 
components. We set Tc(J, K) = 0 everywhere due to the absence of evidence for continuum 
emission or a way to measure it reliably. 


The line shape as a function of frequency, v. is primarily determined by the optical 
depth profile: 


T U (J , K) = t 0 ( J, K) ^2 a i( J > K ) ex P 

i =1 


^—4 In 2 


(v-v c {J,K)-6vi{J,K) 

V Au(J,K) 



(B3) 


where To(J,K) is the opacity in the ( J,K) line, A f(J,K) is the number of hyperhne com¬ 
ponents in the (J, AT) line, i is the index for the hyperhne components, a*(J, Jl) is the scale 
factor for the relative intensity contained in each hyperhne component, n 0 (J, K) is the rest 
frequency of the (.7, K) line, z/ c (J, A') is the Doppler shifted cental frequency of the (J, A') 
line, J, K ) is the difference in frequency of the 7th hyperhne component of the (J, K ) line 


and the main hyperhne component, and Av(J,K) is the FWHM of the (J, K) line (Friesen 


et al. 2009). 


Treating NH 3 as a two-state system, the ratio between the optical depths of two different 
hne structure lines is given by Ho & Townes (1983) as 


t 0 (J',k') = f vq(J',k ') \ 2 / AuQAA^ y 1 Azyjyr) y 1 ( 
r 0 {J,K) \u 0 (J,K)J \Au(J,K) ) \ T ex (J, K) ) [ ^(J, K )\ 2 J \ g(J, K) ) 

( AE(J', K'\ J, K) \ 
6XP l k B T R (J',K'-,J,K))' 

(B4) 


where AE(J', K'\ J , K ) is the diherence in energy of the lower levels of the transitions, 
T R (J',/l / ; J, K ) is the rotation temperature of the lines, |/i(J, 7l)| 2 = H 2 D K 2 / (</(</ + 1)) for 
the electric dipole moment of the molecule //£>, and g(J,K) is the statistical weight given by 


Osorio et al. (2009) as 


g(J, K) 


4(2 J + 1) for K ± 3 or K = 0 
8(2J + 1) for K = 3 and K ^ 0, 


(B5) 


where 3 denotes being a multiple of 3. This accounts for the rotational degeneracy (g u = 
2 J u + 1), the K degeneracy (gx = 1 for K = 0 or gx = 2 for K ^ 0 in symmetric 
top molecules like NH 3 ), and the nuclear spin degeneracy (^nuclear = 4 for K = 3 (ortho- 
NH 3 ) or (/nuclear = 2 for K ^ 3 (para-NH 3 )). We further assume the excitation tem¬ 
perature, the central velocity, and the velocity FWHM are constant for all (J, K ) (thus 
Av{J',K')/Av{J,K) = v Q (J' ,K')/v q (J,K)). 
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From the fit parameters, it is straightforward to calculate the kinetic temperature of 
the gas and the total column density of the NH 3 . The kinetic temperature, T K , is given by 


Tafalla et ah (2004) as 


T k = 7r(2, 2 ; 1,1) 


1 - 


Tr( 2, 2; 1,1) 

42 K 


In ( 1 + 1.1 exp 


16 K 


L T r (2,2;1,1)J 


-l 


(B6) 


and the column density of NH 3 in the (1,1) state is given by Friesen et al. (2009) as 

N( 1 , 1 ) = 


8 ? TV$ 9ll + exp (hUp/ [k B T ex ]) 
c 2 A( 1,1) g 2 1 - exp (hv 0 / [k B T ex ]) 


r v dv, 


(B7) 


where the Einstein A coefficient A(l,l) = 1.68 x 10 ' s 1 (Pickett et ah 1998), g\ = g 2 for 
the statistical weights of the upper and lower energy levels of the (1,1) transition, and 

( \/7rTo(l, l)Az/(l,l) 


r„ dv = 




(B8) 


2\ / \n2 

The total NH 3 column density can then be calculated by considering the partition function 
for NH 3 and using the relation 1V(NH 3 ) = N(l,l)Z tot /Z(l,l) (Friesen et ah 2009). The 
terms of the partition function are calculated by 


Z(J, J ) = (2 J + 1) S(J) exp 


-h[BJ(J+l) + (C - B)J 2 } 

Mk(2, 2; 1,1) 


(B9) 


where the rotational constants are B = 298.117 GHz and C = 186.726 GHz (Pickett et ah 


1998), and S(J) is the statistical weight for ortho- or para-NH 3 . For ortho-NH 3 (J a multiple 


of 3), S(J) = 2, and otherwise S(J) = 1 for para-NH 3 , such as (1,1). Z tot is then just the 
sum of the Z(J, J) terms over all J, which quickly converges. In the rest of this paper, we 
drop the (2, 2; 1,1) notation for simplicity. 


C. Projection Effects 


Many of the quantities computed in this paper depend on the inclination of the elongated 
structures with respect to the line of sight, 6. We cannot make corrections for projection 
effects for individual clumps without a method to measure 6. We can discuss how various 
quantities scale with 9 and, assuming a random distribution of inclinations in space, compute 
typical corrections. Following Peretto et ah (2014), we consider 6 = 45° and 6 = 67° (the 
average inclination angle for randomly sampled filaments is 9 = 67°). The projection effects 
will be quite large for clumps elongated nearly along the line of sight, but we expect that 
only a small fraction of clumps have this property. 
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Masses, column densities, and parameters from extinction mapping and spectral line 
fitting are all relatively insensitive to projection effects. The viewing angle most directly 
affects the measured aspect ratio of the clumps. The “true” aspect ratio, A, is related to the 
aspect ratio we measure by A = Ao(sin#) -1 . This will also manifest itself in the calculation 
of R e ff. If we treat elongated clumps as prolate spheroids that are longer by a factor A 0 
along one axis compared to the other two, then the projected area on the sky also scales as 
(sin#) -1 . The effective radius then scales as R c g oc (sin#) -1 / 2 . We list the parameters that 
depend on the inclination angle and their scaling in Table [7| as well as the correction factor 
for # = 45° and # = 67°. One can see the parameter that is most sensitive to viewing angle 
is £ff, C yi) which is changed by less than a factor of 2 even for # = 45°. 
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Tabic 7. Projection Effects 


Parameter Scaling Correction Factor 

9 = 45° 9 = 67° 


^4o 

oc 

(sin#) -1 

1.41 

1.09 

-Reff 

oc 

(sin #)~ 1,/2 

1.19 

1.04 

^ff,sph 

oc 

(sin 9 )~ 3 ^ 4 

1.30 

1.06 

^ff,cyl 

oc 

(sin#)~ 7//1 

1.83 

1.16 

^vir,sph 

oc 

(sin 9 )~ 1//2 

1.19 

1.04 

^vir,cyl 
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(sin#) -1 

1.41 

1.09 

B C r 

oc 

sin# 

0.71 
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